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ABSTRACT 

We present an imaging study of a sample of eight asymptotic giant branch (AGB) stars in the 
Hi 21-cm line. Using observations from the Very Large Array, we have unambiguously detected Hi 
emission associated with the extended circumstellar envelopes of six of the targets. The detected 
Hi masses range from A/hi ~ 0.015 — 0.055Af Q . The Hi morphologies and kinematics are diverse, 
but in all cases appear to be significantly influenced by the interaction between the circumstellar 
envelope and the surrounding medium. Four stars (RX Lep, Y UMa, Y CVn, and V1942 Sgr) are 
surrounded by detached Hi shells ranging from 0.36 to 0.76 pc across. We interpret these shells as 
resulting from material entrained in a stellar outflow being abruptly slowed at a termination shock 
where it meets the local medium. RX Lep and TX Psc, two stars with moderately high space velocities 
(Kipace > 56 km s _1 ), exhibit extended gaseous wakes (~0.3 and 0.6 pc in the plane of the sky), trailing 
their motion through space. The other detected star, R Peg, displays a peculiar "horseshoe-shaped" 
Hi morphology with emission extended on scales up to ~1.7 pc; in this case, the circumstellar debris 
may have been distorted by transverse flows in the local interstellar medium. We briefly discuss our 
new results in the context of the entire sample of evolved stars that has been imaged in H I to date. 
Subject headings: stars: AGB and post- AGB - stars: winds, outflows - circumstellar matter - radio 
lines: stars 



1. INTRODUCTION 

Stars on the asymptotic giant branch (AGB) undergo 
substantial mass-loss through cool, low-velocity winds 
(V ut ~10 km s _1 ). These winds are one of the domi- 
nant means through which chemically enriched material 
is recycled back into the Galaxy to sustain new gener- 
ations of star and planet formation (e.g., Schroder & 
Scdlmayr 2001; Cristallo et al. 2009; Leitner & Kravtsov 
2011). The distribution of circumstellar debris on the 
largest scales may also be important for the subsequent 
evolution of novae, planetary nebulae and some Type la 
supernovae (e.g., Wang et al. 2004; Mohamed & Pod- 
siadlowski 2007; Wareing et al. 2007a; Chiotellis et al. 
2012; Roy et al. 2012). In addition, the outflows and 
space motions of these mass-losing stars play a role in 
shaping the structure of the interstellar medium (ISM) 
on parsec and sub-parsec scales (Villaver et al. 2002; 
Wareing et al. 2007b; Begum et al. 2010). 

It has long been recognized that the interstellar en- 
vironment of evolved stars can have a profound impact 
on the observed mass, size, shape, and chemical compo- 
sition of their circumstellar envelopes (hereafter CSEs; 
e.g., Smith 1976; Isaacman 1979; Young et al. 1993b; 
Villaver et al. 2002, 2012; Wareing et al. 2006; 2007a, b, 
c; Wareing 2012; see also the review by Stencel 2009). 
Nonetheless, the interfaces through which AGB stars 
interact with their environments have remained poorly 
studied observationally. One problem is that CSEs can 
be enormously extended ( >1 pc), and their chemical 
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compositions change as a function of distance from the 
star as densities drop, and molecules become dissoci- 
ated by the interstellar radiation field or other radiation 
sources. The result is that some of the most frequently 
used observational tracers of CSEs (e.g., CO; SiO, H 2 0, 
and OH masers) do not probe the outermost CSE or its 
interaction zone with the ISM. In some cases, more ex- 
tended material can be traced via imaging observations 
in the far-infrared (FIR; e.g., Young et al. 1993a, b; Zi- 
jlstra & Weinberger 2002; Ueta et al. 2006, 2010; Cox 
et al. 2012a) or in the far-ultraviolet (FUV; Martin et 
al. 2007; Sahai & Chronopoulos 2010), but such data do 
not supply any direct kinematic information. 

Recent studies have demonstrated that the Hi 21- 
cm line of neutral hydrogen is a powerful tool for the 
study of extended CSEs and their interface with their 
larger-scale environments (e.g., Gerard & Le Bertre 2006; 
Matthews & Reid 2007; Libert et al. 2007, 2008, 2010a, 
b; Matthews et al. 2008, 2011a, b, 2012; Gerard et 
al. 2011a, b). Not only is hydrogen the most plenti- 
ful constituent of CSEs, but Hi line measurements pro- 
vide valuable kinematic information and can be directly 
translated into measurements of atomic hydrogen mass. 
Furthermore, because Hi is not destroyed by the inter- 
stellar radiation field, H I measurements are particularly 
well-suited to probing the outermost reaches of CSEs. 

While the number of evolved stars detected in the H I 
line has grown steadily over the past decade thanks to 
ongoing surveys with the Nangay Radio Telescope (NRT; 
Le Bertre & Gerard 2001, 2004; Gerard & Le Bertre 2003; 
Gerard & Le Bertre 2006; Gardan et al. 2006; Libert et 
al. 2008, 2010a; Gerard et al. 2011a, b), the number of 
stars that has been imaged in H I with aperture synthesis 
arrays has remained small (Bowers & Knapp 1987, 1988; 
Matthews & Reid 2007; Matthews et al. 2008, 2011b, 
2012; Le Bertre et al. 2012). Such imaging observa- 
tions offer a crucial complement to single-dish measure- 
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ments, since the high angular resolution of interferome- 
ters provides detailed information on the structure and 
spatially resolved kinematics of the CSE and enables de- 
tailed comparisons with observations of the CSE at other 
wavelengths. In addition, interferometers act as effective 
filters against large-scale Galactic contamination along 
the line-of-sight, allowing stellar signals to be more read- 
ily disentangled from background radiation. 

In this paper we present new H I imaging observations 
of eight mass-losing AGB stars, effectively doubling the 
existing sample. We describe in detail the results for 
each individual star, followed by a discussion of emerging 
trends from the combined sample of evolved stars that 
has been imaged in Hi to date. 

2. SAMPLE SELECTION 

The targets for the present study of circumstellar H I 
comprise a set of eight Galactic AGB stars spanning 
a range in properties, including effective temperature, 
mass-loss rate, variability class, and chemistry. Some ba- 
sic parameters of the sample are summarized in Table 1 . 
None of the stars in our current program is a known or 
suspected binary, although two binary AGB stars, R Aqr 
and Mira, have been observed in H I as part of our earlier 
studies (Matthews & Reid 2007; Matthews et al. 2008). 

Six of our targets (RX Lep, Y UMa, R Peg, Y CVn, 
V1942 Sgr, and TX Psc) were selected from the growing 
sample of AGB stars that has been detected in H I using 
single-dish observations with the NRT (see above). As 
we were interested in imaging examples of both carbon- 
rich and oxygen-rich stars, we selected three from each 
chemistry type, all of which are relatively nearby and 
exhibited relatively strong peak H I signals based on the 
previous NRT measurements ( >0.1 Jy). In the case of 
Y UMa, R Peg, Y CVn and VT942 Sgr, an additional 
motivation for selection was that line-of-sight confusion 
from Galactic emission near the stellar systemic velocity 
is low. Imposing this selection criterion results in the 
selection of stars at high Galactic latitudes that are more 
likely to be Population II stars. In contrast, RX Lep and 
TX Psc exhibit moderate and moderately strong line-of- 
sight confusion, respectively, but were included in our 
sample owing to evidence for elongated H I morphologies 
based on single-dish measurements (Gerard & Le Bertre 
2006). Finally, we downloaded previously unpublished 
data for two stars (IK Tau and AFGL 3099) from the 
Very Large Array (VLA) archive (see Table 2). Both 
of these stars were previously targeted in single-dish H I 
surveys but neither was detected (Zuckerman et al. 1980; 
Gerard & Le Bertre 2006). 

3. OBSERVATIONS AND DATA REDUCTION 

The data that we present here are a combination of 
new and archival Hi 21-cm line observations obtained 
with the VLA of the National Radio Astronomy Obser- 
vatory (NRAO)E All of the data were obtained using the 
most compact (D) configuration (0.035-1.0 km baselines) 
of the VLA, providing sensitivity to emission on scales 
of up to — 15'. The primary beam of the VLA at our ob- 
serving frequency is ~ 31'. The data were reduced and 

5 The National Radio Astronomy Observatory is operated by 
Associated Universities, Inc., under cooperative agreement with 
the National Science Foundation. 



imaged using the Astronomical Image Processing System 
(AIPS). Further details are described below. 

3.1. New VLA Observations 

Our new VLA observations for six stars were obtained 
throughout the transition period during which the orig- 
inal VLA antennas were being retrofitted as part of the 
Expanded Very Large Array (EVLA) upgrade. Through- 
out this period, the array contained fewer than the nom- 
inal 27 antennas (Table 2). As described below, the use 
of this "hybrid" array necessitated special care during 
the calibration process. 

For the observations of R Peg and Y UMa, the VLA 
correlator was used in dual polarization (2AC) mode with 
a 0.78 MHz bandpass, yielding 256 spectral channels with 
3.05 kHz spacing, providing an effective velocity resolu- 
tion of -0.64 km s" 1 . For V1942 Sgr, RX Lep, TX Psc, 
and Y CVn, the VLA correlator was used in 1A mode 
with a 1.5 MHz bandwidth and 512 spectral channels. 
This mode provided a single right circular (RR) polar- 
ization and a 3.05 kHz channel spacing. 

The band center velocities adopted for each star are 
summarized in Table 2. In some cases, the band center 
was offset from the systemic velocity of the star to shift 
strong Galactic emission away from the band edges. To- 
tal integration times spent on each target are also sum- 
marized in Table 2, along with additional details related 
to the observations. 

The observations of the target stars were interspersed 
with observations of a complex gain calibrator approxi- 
mately every 20-25 minutes. Either 3C48 (0137+331) or 
3C286 (1331+305) was used as an absolute flux calibra- 
tor, and in all cases, an additional strong point source 
was observed as a bandpass calibrator (see Table 3). To 
insure that the absolute flux scale and bandpass calibra- 
tion were not corrupted by Galactic emission in the band, 
the flux and bandpass calibrators were each observed 
twice, with the band shifted by sufficient positive and 
negative frequency offsets so as to shift Galactic H I emis- 
sion out of the band. H I survey spectra from Kalberla 
et al. (2005) were used to determine the required fre- 
quency offsets, and the adopted center frequencies are 
summarized in Table 3. Each of the gain calibrators was 
also observed once at each of these offset frequencies to 
insure accurate bootstrapping of the absolute flux scale. 
We estimate that the resulting flux density scale for each 
target star has an uncertainty of <10% (see also H3.3I) . 

3.1.1. Calibration of the R Peg and Y UMa Observations 

Traditionally, the calibration of VLA data sets employs 
a "channel 0" data set formed by taking a vector average 
of the inner 75% of the observing band. However, the 
mismatch in bandpass response between the VLA and 
EVLA antennas causes closure errors in the channel 
data computed this wayQ Furthermore, the hardware 
used to convert the digital signals from the EVLA an- 
tennas into analog signals for the VLA correlator caused 
aliased power in the bottom 0.5 MHz of the baseband 
(i.e., at low- numbered channels) Q This aliasing is of con- 
sequence because of the narrow bandwidth used for our 

e http:/ /www. via. nrao.edu/astro/guides/evlareturn/postproc/ 
7 http: / / www.vla.nrao.edu/astro/guides / evlareturn / aliasing/ 
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TABLE 1 
Properties of the Target Stars 



Name 


o(J2000.0) 


<5(J2000.0) 


I 


b 


K.lsr d 




Spec. 


Var. 


M 










(deg) 


(deg) 


(km s — 1 ) (pc) 


(K) 


Type 


Class 


(M Q yr- 1 ) 


(km s _1 ) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) (7) 


(8) 


(9) 


(10) 


(11) 


(12) 




Oxygen-rich stars 


IK Tau 


03 53 28.87 


+ 11 24 21.7 


177.95 


-31.41 


34.6 245 


2000-3000 


M6e-M10e 


Mira 


4.2 X lO"" 


22.0 


RX Lep 


05 11 22.87 


-11 50 56.72 


212.54 


-27.51 


28.9 149 


3300 


M7 


SRb 


2.0 x 10~ 7 


4.1 


Y UMa 


12 40 21.28 


+55 50 47.61 


126.16 


+61.21 


19.0 385 


3100 


M7II-III 


SRb 


2.6xl0" 7 


5.4 


R Peg 


23 06 39.17 


+10 32 36.10 


85.41 


-44.56 


24.0 400 


2300-2900 


M7e 


Mira 


5.3 x 10" 7 


4.3 


Carbon-rich stars 


Y CVn 


12 45 07.82 


+45 26 24.92 


126.44 


+71.64 


21.1 272 


2730 


C5.4J 


SRb 


1.7 x 10-'' 


7.8 


V1942 Sgr 


19 19 09.59 


-15 54 30.03 


21.61 


-13.19 


-33.7 535 


2960 


C 


Lb 


(6.1, 1.0) x 10~ 7 


17.5,5.0 


AFGL 3099 


23 28 17.50 


+ 10 54 36.9 


92.23 


-46.94 


46.6 1500 


1800-2000 


C 


Mira 


1.3 x 10~ 5 


10.1 


TX Psc 


23 46 23.51 


+03 29 12.52 


93.28 


-55.59 


12.2 275 


3115 


C5II 


Lb 


3.1 x 10" 7 


10.7 



NOTE. — Units of right ascension are hours, minutes, and seconds. Units of declination are degrees, arcminutcs, and arcseconds. Explanation of columns: (1) star name; (2) & (3) 
right ascension and declination (.12000.0); (4) & (5) Galactic coordinates; (6) systemic velocity relative to the Local Standard of Rest (LSR); (7) adopted distance in parsccs; (8) stellar 
effective temperature in kclvins; (9) spectral type; (10) variability class; (11) mass-loss rate in solar masses per year; two values arc quoted in cases of multi-component molecular 
line profiles (12) outflow velocity in km s ~ 1 as determined from CO observations. Coord inates an d spectral classifications were taken from SIMBAD I http : //simbad . harvard .cdu I. 
References and additional descriptions for the other parameters are provided throughout E|4. If H4.8| 



TABLE 2 

Summary of VLA Observations 



Star 


Obs. date 


No. 


Mode 


No. 


Bandwidth 


v u v 2 


Vccnt 


Au 


Ad 


t 






antennas 




pol. 


(MHz) 


(km s" 1 ) 


(km s —1 ) 


(kHz) 


(km s" 1 ) 


(hours 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


IK Tau* 


1989-Dec-18 


27 


2AC 


2 


0.39 


-4.6,72.6 


34.0 


12.2 


2.6 


3.97 


RX Lep 


2009-Dec-30/31 


22 


1A 


1 


1.5 


-135.3,193.2 


29.0 


3.05 


1.29 


3.85 




2010-Jan-l/2 


22 


1A 


1 


1.5 


-135.3,193.2 


29.0 


3.05 


1.29 


2.75 


Y UMa 


2007-Apr-27 


26 


2AC 


2 


0.78 


-62.8.100.8 


19.0 


3.05 


0.64 


3.35 




2007-Apr-28 


26 


2AC 


2 


0.78 


-62.8.100.8 


19.0 


3.05 


0.64 


3.45 


R Peg 


2007-Apr-8 


26 


2AC 


2 


0.78 


-81.8,81.8 


0.0 


3.05 


0.64 


3.40 




2007-May-4 


26 


2AC 


2 


0.78 


-81.8,81.8 


0.0 


3.05 


0.64 


3.42 


Y CVn 


2010-Jan-l/2 


22 


1A 


1 


1.5 


-143.1,183.5 


21.1 


3.05 


1.29 


5.10 




2010-Jan-2/3 


21 


1A 


1 


1.5 


-143.1,183.5 


21.1 


3.05 


1.29 


0.90 


V1942 Sgr 


2009-Dec-5 


23 


1A 


1 


1.5 


-196.2,132.2 


-32.0 


3.05 


1.29 


2.93 




2009-Dec-7 


22 


1A 


1 


1.5 


-196.2,132.2 


-32.0 


3.05 


1.29 


2.88 


AFGL 3099* 


1989-Dec-17 


27 


2AC 


2 


0.19 


26.7,65.3 


46.0 


6.10 


1.29 


4.23 


TX Psc 


2009-Dec-5/6 


23 


1A 


1 


1.5 


-151.3,177.3 


13.0 


3.05 


1.29 


0.80 




2009-Dec-7/8 


22 


1A 


1 


1.5 


-151.3,177.3 


13.0 


3.05 


1.29 


0.82 




2010-Jan-4/5 


21 


1A 


1 


1.5 


-151.3,177.3 


13.0 


3.05 


1.29 


4.00 



NOTE. — Explanation of columns: (1) star name; (2) observing date; (3) number of available antennas; (4) VLA correlator mode; (5) number of circular 
polarizations; (6) total bandwidth; (7) minimum and maximum LSR velocity covered by the observing band; (8) LSR velocity at band center; (9) channel width; 

(10) effective velocity resolution of the data presented in this work; all data with the exceptions of Y UMa and R Peg were Hanning smoothed prior to analysis; 

(11) total on-sourcc integration time. 

* . . 

Archival data (sec iTOl . 



observations. However, it affects EVLA-EVLA baselines 
only, as it does not correlate on ELVA-VLA baselines. 

To mitigate the above effects, we used the following 
modified approach to the gain calibration for the R Peg 
and Y UMa data sets. After applying the latest avail- 
able corrections to the antenna positions and perform- 
ing an initial excision of corrupted data, we computed 
and applied an initial bandpass calibration to our spec- 
tral line data to remove closure errors on VLA-EVLA 
baselines. EVLA-EVLA baselines were flagged during 
this step and unflagged afterward. After applying the 
bandpass solutions, we then computed a new frequency- 
averaged (channel 0) data set (using the inner 75% of the 
band) for use in calibrating the frequency-independent 
complex gains. After flagging the EVLA-EVLA base- 
lines on the channel data, we solved for the frequency- 
independent portion of the complex gain and transferred 
these solutions to the full spectral line data set. Lastly, 
we applied time-dependent frequency shifts to the data 



to compensate for changes caused by the Earth's motion 
during the course of the observations. 

3.1.2. Calibration of the RX Lep, Y CVn, V1942 Sgr, and 
TX Psc Observations 

During the observations of RX Lep, Y CVn, V1942 Sgr, 
and TX Psc, there were too few non-EVLA antennas 
left in the array to utilize the calibration strategy de- 
scribed in §3.1.11 Instead, we increased our observing 
bandwidth by a factor of two (at the cost of limiting 
ourselves to a single polarization product) to provide a 
suitable amount of bandwidth uncorrupted by aliasing. 
We then employed the following calibration approach. 

After applying the latest available antenna position 
corrections and performing an initial excision of cor- 
rupted data, we computed and applied a bandpass cal- 
ibration to the spectral line data to remove closure er- 
rors on VLA-EVLA baselines. All baselines were used in 
the calculation, but the bandpass was normalized using 
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TABLE 3 

Calibration Sources 



Source 


a(J2000.0) 


<5(J2000.0) 


Flux Density (Jy) 


u (MHz) 


Date 




ui o / 4i.zyy4 


iqq no Qc; i on 
~\-oo Uy oo.LoZ 


lO.ooOl 


1 A 1 Q K 


onnv A r*r Q 

ZUU / -Apr-o 








10. oO I 


1 A on e; 
14ZU.D 


onnv A T^t- q 
zUU / -i\pi-o 








10.oo4 f 


1 /1 1 n c; 


onn7 l\/r^., r a 
zUU /-IVIay-4 








1 ^ * 

lU.OU ( 1 


1420 5 


zuu / -iviay-^i 








1 5 8905* 


1418.9 


2009-Dec-6 








1 5 8691 * 


1421.3 


2009-Dec-6 








1 5 8905* 


1418.9 


2009-Dec-8 








15.8691* 


1421.3 


2009-Dec-8 








15.8885* 


1419.1 


2009-Dec-31 








15.8701* 


1421.5 


2009-Dec-31 








15.8885* 


1419.1 


2010-Jan-2 








15.8701* 


1421.5 


2010-Jan-2 








15.8868* 


1419.3 


2010-Jan-5 








15.8696* 


1421.2 


2010-Jan-5 


0022+002 b 


00 22 25.42 


+00 14 56.16 


2.92+0.06 


1420.2 


2009-Dec-6 








3 02+0 07 


1420.2 


2009-Dec-8 








2 73+0 05 


1420.2 


2010- Jan-2 


0409+122 c 


04 09 22.00 


+12 17 39.84 


0.814+0.002 


1420.7 


1989-Dec-18 


04^Q4-094 d 


04 ^9 n 1 ^ 


4-09 9Q "il 17 


i 794-0 04. 


1420.1 










i 7Q4-0 OS 


1420.1 


2010-Jan-2 


1227+365 6 


12 27 58.72 


+36 35 11.82 


1.97+0.04 


1420.4 


2010-Jan-2 








2 07+0 06 


1420.4 


2010-Jan-3 


1 OE;o_l_c:fie; f 
IZOZ-hODO*' 


1 9 PiO Oft 3fi 


+ 00 o4 iy.4i 


o onin no 


i a on 3 


onn7 A rxv o? 
zuu / -/\pr-z / 










1 a on Q 


0n07 A rxr- OS 

zuu / -/\pr-zo 




1 Q Q1 AC 9S7Q 




1 fi o^qzl* 


1420 3 


1 QSQ T)or 1 & 
J- yoy-iJct,-io 








14.7317* 


1419.3 


2007-Apr-27 








14.7218* 


1421.3 


2007-Apr-27 








14.7317* 


1419.3 


2007-Apr-28 








14.7218* 


1421.3 


2007-Apr-28 








14.7305* 


1419.5 


2010-Jan-2 








14.7193* 


1421.8 


2010-Jan-2 








14.7305* 


1419.5 


2010-Jan-3 








14.7193* 


1421.8 


2010-Jan-3 


1411+522 h 


14 11 20.64 


+52 12 09.14 


21.63+0.08 


1419.3 


2007-Apr-27 








21.72+0.06 


1421.3 


2007-Apr-27 








21.84+0.05 


1419.3 


2007-Apr-28 








21.78+0.05 


1421.3 


2007-Apr-28 








21.83+0.20 


1420.4 


2010-Jan-2 








22.01+0.73 


1420.4 


2010-Jan-3 


1911-201 1 


19 11 09.65 


-20 06 55.10 


2.17+0.04 


1420.5 


2009-Dec-5 








2.24+0.08 


1420.5 


2009-Dec-7 


2253+16P 


22 53 57.74 


+16 08 53.56 


14.71+0.04 


1419.5 


2007-Apr-8 








14.69+0.04 


1421.5 


2007-Apr-8 








14.46+0.05 


1419.5 


2007-May-4 








14.51+0.04 


1421.5 


2007-May-4 








14.13+0.10 


1418.9 


2009-Dec-7 








13.42+0.17 


1419.5 


2010-Jan-4 








13.31+0.14 


1421.8 


2010-Jan-4 


2255+132 fc 


22 55 03.73 


+13 13 33.01 


2.65+0.01 


1420.5 


2007-Apr-8 








2.66+0.02 


1420.5 


2007-May-4 


2330+110' 


23 30 40.85 


+11 00 18.70 


1.117+0.003 


1420.0 


1989-Dec-17 



NOTE. — Units of right ascension are hours, minutes, and seconds, and units of declination arc degre es, a rcminutcs, and arcseconds. i>q 
is the frequency at which the flux density in the fourth column was computed. For the archival data f H3 . 2\ , the phase center coordinates 
and phase calibrator names have been translated from B1950.0 to J2000.0 coordinates. 
* 

Adopted flux density, computed according to the VLA Calibration Manual (Perley & Taylor 2003). 
a Primary flux calibrator and bandpass calibrator for RX Lcp, R Peg, V1942 Sgr, and TX Psc 
k Secondary gain calibrator for TX Psc 
C Secondary gain calibrator for IK Tau 

Secondary gain calibrator for R.X Lcp 
C Secondary gain calibrator for Y CVn 
f Secondary gain calibrator for Y UMa 

S Primary flux calibrator for IK Tau, Y UMa, Y CVn, and AFGL 3099: bandpass calibrator for Y UMa and Y CVn 
h Bandpass calibrator for Y UMa and Y CVn 
1 Secondary gain calibrator for V1942 Sgr 
J Bandpass calibrator for R Peg and TX Psc 

Secondary gain calibrator for R. Peg 
1 Secondary gain calibrator for AFGL 3099 
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TABLE 4 

Deconvolved Image Characteristics 



Source 


Type 


n 


Taper 


#FWHM 


PA 


rms 


Continuum velocities 


Clean Boxes? 








(kA,kA) 


("x") 


(degrees) 


(mjy beam -1 ) 


(km s- 1 ) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


IK Tau 


line 


+5 




69'.'2x51'.'7 


32.1 


1.4 


57.2 to 72.6 


no 


IK Tau 


cont. 


+ 1 




69'.'2x51'.'7 


32.1 


0.61 


57.2 to 72.6 


yes 


RX Lep 


line 


+5 




76'.'2x47'.'9 


10.6 


2.0 


-37.3 to -121.0; 34.8 to 168.7 


yes 


RX Lep 


cont. 


+ 1 




67'.'5x43'.'8 


10.7 


0.4 


-37.3 to -121.0; 34.8 to 168.7 


yes 


Y UMa° 


line 


+1 




61'.'2x52'.'l 


51.2 


1.7 


-41.5 to -15.1; 0.96 to 12.5; 24.2 to 88.6 


yes 


Y UMa° 


line 


+5 




72'.'lx58'. / l 


54.9 


1.7 


-41.5 to -15.1; 0.96 to 12.5; 24.2 to 88.6 


yes 


Y UMa 


cont. 


+1 




61'.'7x52'.'4 


51.4 


0.18 


-41.5 to -15.1; 0.96 to 12.5; 24.2 to 88.6 


yes 


R Peg" 


line 


+5 




62'.'9x56'.'9 


60.9 


1.7 


-68.9 to -32.8; 35.4 to 57.3 


no 


R Peg 


cont. 


+1 




56'.'9x67'.'2 


67.2 


0.24 


-68.9 to -32.8; 35.4 to 57.3 


yes 


Y CVn 


line 


+1 




53'.'6x44'.'4 


30.1 


2.2 


-125.1 to -25.9; 57.8 to 100.3; 124.8 to 150.5 


yes 


Y CVn 


line 


+5 




60'.'9x48'.'6 


32.8 


2.1 


-125.1 to -25.9; 57.8 to 100.3; 124.8 to 150.5 


yes 


Y CVn 


line 


+5 


2,2 


106'.'9x83'.'5 


38.6 


2.6 


-125.1 to -25.9; 57.8 to 100.3; 124.8 to 150.5 


yes 


Y CVn 


cont. 


+1 




53'.'6x44'.'4 


30.0 


0.48 


-125.1 to -25.9; 57.8 to 100.3; 124.8 to 150.5 


yes 


V1942 Sgr 


line 


+5 




89'.'5x49'.'4 


12.8 


2.7 


-175.6 to -50.6 


yes 


V1942 Sgr 


cont. 


+ 1 




79'.'6x45'.'7 


13.8 


0.48 


-175.6 to -50.6 


yes 


AFGL 3099 


line 


+5 




63'.'6x51'.'9 


15.6 


1.7 


26.7 to 39.5; 52.4 to 65.3 


no 


AFGL 3099 


cont. 


+ 1 




56'.'4x47'.'9 


13.8 


0.45 


26.7 to 39.5; 52.4 to 65.3 


yes 


TX Psc 


line 


+5 




71'.'6x48'.'8 


5.0 


2.0 


-130.6 to -25.0; 75.5 to 150.2 


no 


TX Psc 


line 


+5 


2.2,2.2 


104'.'2x78'.'6 


4.2 


2.2 


-130.6 to -25.0; 75.5 to 150.2 


no 


TX Psc 


cont. 


+ 1 




63'.'6x45'.'6 


6.9 


0.29 


-130.6 to -25.0; 75.5 to 150.2 


yes 



NOTE. — Explanation of columns: (1) target name; (2) indication of whether the image contains line or continuum emission: the continuum images comprise a single spectral 
channel representing an average of the line-free portions of the band: (3) AIPS robust parameter used in image dcconvolution; 7?. — + 5 is comparable to natural weighting: (4) 
Gaussian taper applied in u and v directions, expressed as distance to 30% point of Gaussian in units of kilolambda; (5) FWHM dimensions of synthesized beam; (6) position angle 
of synthesized beam (measured east from north); (7) rms noise per channel (Ict); (8) LSR velocities of spectral channels used for continuum subtraction {line data) or that were 
averaged to compute a continuum image; (9) indication of whether or not clean boxes were used during image deconvolution . 

a Data were not Hanning smoothed for the present analysis. 
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only channels 163 through 448, thus excluding the por- 
tion of the band affected by aliasing. We next computed 
a frequency-averaged channel data set for use in cali- 
brating the frequency-independent complex gains, using 
a vector average of channels 163 through 448. Following 
gain calibration, we applied time-dependent frequency 
shifts to the data to compensate for changes caused by 
the Earth's motion during the course of the observa- 
tions. Lastly, we applied Hanning smoothing in veloc- 
ity and discarded every other channel, yielding a 256 
channel data set with an effective velocity resolution of 
-1.3 km s _1 . 

Issues affecting two stars merit special comment. The 
data for RX Lep were heavily contaminated by narrow- 
band radio frequency interference (RFI), necessitating 
considerable channel-dependent flagging. This resulted 
in a loss of continuum sensitivity, although fortunately, 
little RFI was present near the systemic velocity of the 
star. In the case of V1942 Sgr, a portion of the data 
taken on 2009 December 7 between 23:07:25 UT and 
23:25:46 UT were affected by a bug in the VLA archive 
software that resulted in the frequencies assigned to each 
of the spectral channels being mirror reversed relative to 
their respective channel numbers. The affected data were 
corrected using the AIPS task FLOPM. 

3.1.3. Continuum Subtraction 

Prior to imaging the line data for all of the newly ob- 
served target stars, we performed a continuum subtrac- 
tion using a linear fit to the real and imaginary compo- 
nents of the visibilities using the AIPS task UVLIN. The 
portions of the band that were determined to be line-free 
channels and were used for these fits are summarized in 
Column 7 of Table 4. Although the spectral shape of 
the aliased portion of the continuum as measured toward 
our continuum calibrators was better approximated by a 
higher order polynomial, there was generally insufficient 
continuum signal in our target fields to adequately con- 
strain higher order fits. 

3.2. Archival VLA Data 

H I data for two stars (AFGL 3099 and IK Tau) were 
retrieved from the VLA archive. These data were origi- 
nally taken as part of program AK237 in late 1989 (Prin- 
cipal Investigator: G. Knapp) and to our knowledge are 
previously unpublished. 

The data for IK Tau comprise 32 frequency channels 
with a channel spacing of 12.2 kHz (2.6 km s" 1 ) after on- 
line Hanning smoothing, resulting in a total bandwidth 
of 0.39 MHz (—82 km s _1 ) in dual circular polarizations. 
The AFGL 3099 data comprise 32 frequency channels 
with a channel spacing of 6.1 kHz (1.3 km s _1 ) after on- 
line Hanning smoothing, providing a total bandwidth of 
0.195 MHz (—41 km s _1 ), also with dual circular polar- 
ization. Additional details are provided in Table 2. 

We processed both of the archival data sets in a similar 
manner. After flagging corrupted data, we used a chan- 
nel data set to calibrate the antenna-based, frequency 
independent complex gains. The absolute flux level 
was established using observations of 3C48 (0137+331). 
We adopted the VLA-determined flux density values for 
3C48 from the year 1990, as calculated by the AIPS task 
SETJY (see Table 3). 



Because the 3C48 observations suffered from Hi ab- 
sorption across much of the band, we estimate that the 
uncertainty in the absolute flux scale for these data may 
be as high as —20%. This Hi absorption also rendered 
3C48 unsuitable for use as a bandpass calibrator, hence 
no bandpass calibration was applied to the archival data. 
This had minimal impact on the data quality given the 
narrow bandwidths and the fact that the channels at the 
band edges had already been discarded by the VLA on- 
line system. 

To create a continuum-subtracted spectral data set for 
each of the archival sources, we produced a clean image of 
the continuum using an average of the line-free channels 
(see Table 4) and then subtracted the clean components 
from the line data in the visibility domain using the AIPS 
task UVSUB. 

3.3. Imaging the VLA Data 
3.3.1. Continuum Data 

We produced an image of the 21-cm continuum emis- 
sion in the field of each target star using the line-free 
portion of the band (see Table 4). No continuum emis- 
sion was detected from any of the stars. For the six 
newly observed target fields f ^3.1|) . we compared the 
measured flux densities of several continuum sources in 
each field with values measured from NRAO VLA Sky 
Survey (NVSS) images (Condon et al. 1998). In all cases 
we found agreement to within better than 10%. For each 
target, we also searched for Hi absorption against the 
brightest continuum sources in the field surrounding each 
star, but in no case did we detect statistically significant 
absorption. 

3.3.2. Line Data 

We imaged the VLA Hi line data (both new and 
archival) using the standard AIPS CLEAN deconvolu- 
tion algorithm. For each target, we produced data cubes 
using various weighting schemes. The characteristics of 
those used for the present analysis are summarized in 
Table 4. 

4. RESULTS 

In £14.1114.81 we present the results of our most recent 
Hi imaging study of AGB stars. In six cases we have 
unambiguously detected H I emission associated with the 
circumstellar envelopes of the target stars. For the other 
two targets we derive upper limits on the Hi content 
of their CSEs. We discuss the VLA data for each star 
and their implications for understanding the mass loss 
history of each star in the following subsections. In §5 we 
discuss emerging trends from the entire sample of AGB 
and related stars that has been imaged in H I to date. 

4.1. IK Tauri (IK Tau) 

IK Tauri (IK Tau=NML Tau) is an extremely red, 
large-amplitude Mira variable with an oxygen-rich chem- 
istry and a period of 462 days (Le Bertre 1993). The 
large amplitude of its variability corresponds to varia- 
tions in its effective temperature from —2000 K to 3000 K 
(Dyck et al. 1974; Gerard & Le Bertre 2006). Following 
Le Bertre & Winters (1998), we adopt a distance to the 
star of 245 pc based on the period-luminosity relation of 
Feast (1996). 
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Fig. 1. — Hi channel maps bracketing the systemic velocity of IK Tau, derived from naturally weighted VLA data. The spatial resolution 
is ~ 69" X 52". A star symbol indicates the stellar position. Contour levels are (-8.4, -6, -4.2, -3, 3... 8.4) X 1.4 mjy beam -1 . The lowest 
contour is ~ 3<r and negative contours are shown as dashed lines. The field-of-view shown is comparable to the VLA primary beam. 



IK Tau has been extensively studied at multiple wave- 
lengths and is known to have a rich molecular chem- 
istry (e.g., Milam et al. 2007; Decin et al. 2010; 
Claussen et al. 2011). Based on CO(l-O) observations, 
Knapp & Morris (1985) derived a stellar systemic ve- 
locity of V^lsr = 34.6 ± 0.4 km s _1 , a mass-loss rate 
M = 4.2 x 10~ 6 M Q yr _1 (scaled to our adopted dis- 
tance), and an outflow velocity V ont — 22.0 ±0.6 km s _1 . 

4.1.1. VLA Results for IK Tau 

Hi channel maps bracketing the systemic velocity of 
IK Tau are presented in Figure [0 The velocity range 
shown roughly corresponds to the velocities over which 
CO emission has been previously detected in the CSE. 
We find no Hi emission features with > 3ct significance 
coincident with the stellar position in any of the spectral 
channels with Vlsr > 23.7 km s _1 . For velocities Vlsr < 
28.7 km s _1 , the data suffer from confusion from line-of- 
sight Galactic emission that is poorly spatially sampled 
by the VLA, producing a "mottled" appearance in the 
channel images. There is no evidence that any of the 
emission seen at these velocities is related to the CSE of 
IK Tau. 



In Figure we present a spatially integrated H I spec- 
trum derived by summing each plane of the spectral line 
data cube within a circular aperture of radius 85", cen- 
tered on the star. This aperture was chosen to match the 
extent of the FIR-emitting shell detected by Cox et al. 
(2012a; see H4.1.2p , Consistent with Figure[I] the result- 
ing H I spectrum is contaminated by Galactic emission 
at velocities Vlsr <24 km s _1 . However, near the stellar 
systemic velocity, we see marginal evidence for a weak, 
broad spectral feature that is distinct from the Galactic 
emission component at lower velocities. Using the spec- 
tral data spanning 26.3 < Vlsr < 72.4 km s _1 , we per- 
formed a Gaussian fit to this feature. The best-fit Gaus- 
sian is overplotted on Figure it has a peak amplitude 
of 5.6±1.5 mjy, a velocity centroid of 34.4±2.3 km s _1 , 
a FWHM velocity width of 17.6±2.9 km s~\ and a DC 
offset term of -0.0013±0.0008 Jy. The peak amplitude 
of the Gaussian is roughly 3.5 times the RMS noise in 
the line-free portion of the spectrum, and the velocity 
centroid agrees with the stellar systemic velocity derived 
from CO observations. Despite this, we consider this 
feature to be of marginal significance given the combi- 
nation of several factors. The blend with the Galactic 
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Fig. 2. — Spatially integrated Hi spectrum of IK Tau derived 
from naturally weighted VLA data. The flux density in each chan- 
nel was derive d by in tegrating within a circular aperture with ra- 
dius 85" (see 94.1.1^ . The vertical bar indicates the systemic ve- 
locity of the star derived from CO observations. Spectral channels 
with Vlsr ^26 km s are contaminated by Galactic emission 
along the line-of-sight. The dotted line is a Gaussian fit to the 
data over the velocity range 26.3 < Vlsr ^ 72.4 km s —1 ; we con- 
sider this feature to be of marginal significance (see i|4.1,ljl . 

emission on the blue edge of the profile makes our Gaus- 
sian fit parameters more uncertain than the formal errors 
would indicate, and the emission giving rise to this possi- 
ble spectral feature cannot be unambiguously associated 
with the CSE of IK Tau in any of the individual channel 
images. Finally, our data cube contains too few line- 
free channels to permit a robust characterization of the 
intrinsic instrumental band shape. Deeper H I imaging 
observations with higher spectral resolution and a wider 
band could help to overcome these issues. 

For the present analysis we use the integral of the 
Gaussian fit in Figure [2] to provide a ~ 3(7 up- 
per limit to the integrated Hi flux density from the 
CSE of IK Tau within a 85" (0.10 pc) radius of 
/ Shi dv <0.10 Jy km s _1 . This translates to an Hi 
mass of Mm < 1.4 x 10 _3 Mq at our adopted distance. 

4.1.2. Discussion of IK Tau Results 

IK Tau was previously targeted in single-dish H I sur- 
veys by Zuckerman et al. (1980) and Gerard & Le Bertre 
(2006) but was not detected. IK Tau is a challenging 
target for single-dish studies owing to the strong line-of- 
sight contamination across the velocity range Vlsr ~— 40 
to 40 km s _1 , which encompasses the stellar systemic ve- 
locity (see Gerard & Le Bertre 2006). 

With the VLA, we almost entirely filter out the con- 
taminating emission for Vlsr >31.4 km s _1 , allowing us 
to place more robust constraints on the total amount 
of atomic hydrogen in the CSE. However, we note some 
important caveats to the Hi mass upper limit that we 
derived in the previous section. First, the rather coarse 
spectral resolution employed for the IK Tau observa- 
tions (^2.6 km s" 1 ) may have diminished our sensitiv- 
ity to emission components with narrow velocity widths 
(AV ~ 1 km s _1 ). Such components may be present 
even when the global (spatially integrated) H I line width 
is several times larger. A second caveat is that the gas 
temperature in CSEs is expected to depend strongly on 



the mass-loss rate (Jura et al. 1988; Sahai 1990), with 
cooler CSE temperatures expected for stars with higher 
values of M. Thus a significant fraction of the stellar 
wind may be molecular (Glassgold & Huggins 1983) and 
most of the remaining fraction of the atomic hydrogen 
in the CSE of IK Tau may be simply too cold to de- 
tect in emission. This argument may also be relevant to 
the case of AFGL 3099 ( §4.7j) . another star that remains 
undetected in H I emission despite its high mass loss rate. 

Recent Herschel imaging at 70/im and 160/mi by Cox 
et al. (2012a) revealed that IK Tau is surrounded by ex- 
tended infrared emission with radius ~ 85" that appears 
to be consistent with a bow shock and related structures. 
The total mass of the CSE estimated from these FIR 
measurements is ~ 0.037Mq (assuming a gas-to-dust ra- 
tio of 200; Cox et al. 2012b), roughly 20 times higher 
than the upper limit we derive from our H I measure- 
ments (after scaling by a factor of 1.4 to account for the 
mass of helium) . 

4.2. RX Leporis (RX Lep) 

RX Leporis (RX Lep) is an oxygen-rich AGB star that 
was reclassified by Samus et al. (2004) from an irregular 
variable to a semi-regular variable of type SRb. We adopt 
here a distance of 149 pc based on the Hipparcos parallax 
of 6.71±0.44 mas (van Leeuwen 2007). 

Libert et al. (2008) presented CO(2-l) observations 
for RX Lep from which they derived a stellar systemic 
velocity of V+,lsr = 28.9±0.1 km s _1 , a mass-loss rate 
M 2.0 x 10~ 7 M e yr" 1 (scaled to our adopted dis- 
tance), and an outflow velocity for the wind V ou t — 
4.2 ±0.1 km s _1 . The proper motions measured by Hip- 
parcos are 31.76±0.58 mas yr _1 in right ascension and 
56.93±0.50 mas yr _1 in declination, implying a peculiar 
space motion of V spacc «56.6 km s _1 along a position 
angle of 24°6 (see Table 5). 

RX Lep has an effective temperature of ^3300 K 
(Dumm & Schild 1998), and using the bolometric cor- 
rection from Le Bertre et al. (2001), its luminosity is 
5300 Lq at our adopted distance (Libert et al. 2008). 
The initial mass of RX Lep is believed to be in the range 
2.5-4.0 M Q (Mennessier et al. 2001). As discussed fur- 
ther below, RX Lep was previously detected in Hi by 
Libert et al. (2008) using single-dish observations from 
the NRT. 

Technetium lines have not been detected in RX Lep 
(Little et al. 1987; Lebzelter & Hron 1999), implying 
that either RX Lep is on the early AGB (E-AGB) or 
that it is still near the beginning of the thermal pulsing 
(TP) AGB phase and has not yet undergone dredge-up 
(e.g., Straniero et al. 1997). The mass- loss rate and wind 
outflow speed for the star are also consistent with values 
expected near the beginning of the TP- AGB (Vassiliadis 
& Wood 1993). As described below ( §4.2.2 j) . placement 
of RX Lep near the start of the TP-AGB is further sup- 
ported by our new Hi observations. 

4.2.1. VLA Results for RX Lep 

Hi channel maps for RX Lep derived from our VLA 
data are presented in Figure [3] and an H I total intensity 
map derived from the same data is shown in Figure 2) 
Together these figures reveal the detection of an extended 
H I emission nebula surrounding the position of RX Lep 
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Fig. 3. — Hi channel maps bracketing the systemic velocity of RX Lep, derived from naturally weighted VLA data. The spatial resolution 
is ~ 76" X 48". A star symbol indicates the stellar position. Contour levels are (— 6[absent],— 4.2,— 3,3,4.2,6,8.4,12) X 2.0 mjy beam -1 . The 
lowest contour is ~ 3a. The channels shown correspond to the range of LSR velocities over which CO(2-l) emission was detected by Libert 
et al. 2008. The field-of-view displayed is comparable to the VLA primary beam. 



whose morphology and velocity distribution imply an 
unambiguous association with the star. The H I emis- 
sion is most prevalent in the two velocity channels that 
bracket the stellar systemic velocity (Vlsr = 28.4 km s _1 
and 29.6 km s -1 ), but the emission can be traced from 
25.8 km s -1 to 32.2 km s _1 . 

Additional H I emission is also detected near the posi- 
tion of RX Lep in velocity channels from 16.8 km s _1 
to 24.5 km s" 1 . However, the peak intensity of this 
emission is centered a few arcminutes southeast of the 
star, and it reaches its peak spatially-integrated intensity 
(-1.2 Jy) at Vlsr =20.6 km s" 1 (see Figure 0. This 
is outside the velocity range over which molecular line 
emission has been detected in the CSE of RX Lep. Fur- 
ther, the borders of this emission are not well defined; 
at low surface brightness levels they begin to blend with 
large-scale bands of positive and negative emission that 
stretch across the corresponding channel images. To- 
gether these properties suggest that this emission is not 
related to the CSE of RX Lep, but instead is part of an 
interstellar cloud or cloud complex along the line-of-sight 



that is spatially extended and poorly spatially sampled 
by the VLA. This interpretation is consistent with evi- 
dence for confusion over this velocity range seen in the 
single-dish spectra of Libert et al. (2008). Based on simi- 
lar arguments, the H I emission seen northwest of RX Lep 
(near a J2 ooo = 05 h 10 m 39 s , <5 J2 ooo = -11°39 / 16") is also 
unlikely to be related to the star. 

A defining characteristic of the H I emission associated 
with RX Lep is its "cracked egg" morphology. As seen 
in Figure 21 the velocity-integrated emission exhibits a 
distinctly ovate shape with its semi-major axis aligned 
north-south. The extent of the "egg" is ^6'.7 (0.30 pc) 
east-west by ~9'8 (0.42 pc) north-south. The star itself 
lies ~ 30" north of the geometric center of the egg. The 
H I distribution across the egg appears clumpy, and Fig- 
ure [5] reveals two enhancements in the Hi column density 
at its north and south "poles" . In contrast, a pronounced 
dearth of emission is observed along the egg's east-west 
equator. 

Another striking feature of RX Lep's circumstellar H I 
emission is the presence of a narrow, highly extended (~ 
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Fig. 4. — Hi total intensity map of RX Lep derived from naturally weighted data. The spatial resolution is ~ 76" X 48". The map was 
constructed from emission spanning LSR velocities 25.8 to 32.2 km s — 1 . The red line indicates the direction of space motion of the star. 
Intensity levels are to 70 Jy beam -1 m s _1 and contour levels are (1.8,3.5,3.5,7,10) X 7.3 Jy beam -1 m s — 1 . To minimize noise in the 
map, data at a given point were blanked if they did not exceed a 2.5a threshold after smoothing by a factor of three spatially and spectrally. 



7' or 0.30 pc long) tail of emission that stretches almost 
due south of the star. As seen in Figure 21 the tail in- 
cludes several compact knots of emission. There are also 
hints that the tail may have a kinked or corkscrew shape, 
although this cannot be confirmed from the present data. 

Two versions of the spatially integrated H I spectrum 
of RX Lep are shown in Figure \5\ The first (shown as 
a solid histogram) is derived by integrating the emission 
within each spectral channel within a fixed aperture. Be- 
cause of the presence of confusion blueward of the stellar 
systemic velocity, we have also derived a second spectrum 
(plotted as a dashed histogram) by summing the emission 
in irregularly shaped "blotches" surrounding the stellar 
position in each channel over the velocity range at which 
circumstellar CO emission has been detected. The blotch 
peripheries were defined by 3a contours and encompass 
the emission comprising the tail and shell structures vis- 
ible in Figure 2J 

Based on a two-component Gaussian fit to the "fixed 
aperture" spectrum in Figure [5l we derive an H I line 
centroid for the RX Lep emission of Vlsr = 28.4 ± 



0.1 km s , in good agreement with the value determined 
from CO observations. From our Gaussian fits we also 
find a peak flux density of 0.58±0.02 Jy and a FWHM 
Hi linewidth AV = 5.0±0.3 km s _1 . The latter is larger 
than the value of 3.8±0.1 km s _1 previously reported 
by Libert et al. (2008). From our fit we compute an 
integrated Hi flux density for the CSE of RX Lep of 
/ S v dv = 3.09±0.22 Jy km s _1 . Alternatively, using the 
"blotch" spectrum, we measure J S v dv — 2.77±0.04 Jy 
km s . Because the uncertainty is dominated by the 
systemic uncertainties in determining the intrinsic line 
profile shape in the presence of confusion, we adopt 
as our best estimate a mean of our two measurements: 
/ S v dv — 2.93 ± 0.22 Jy km s _1 . This translates to an 
Hi mass of M H i ~ 0.015M Q . This is roughly 2.5 times 
larger than measured by Libert et al. (2008) from single- 
dish mapping. The reason for this discrepancy is unclear, 
although part of the difference may arise from the man- 
ner in which the RX Lep line emission was deblended 
from the neighboring Galactic signal. Libert et al. used 
polynomial baseline fits, while here we use spatial segre- 
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Fig. 5. — Spatially integrated Hi spectrum of RX Lep. The 
stellar systemic velocity derived from molecular line observations is 
indicated by a vertical bar. The thick histogram shows a spectrum 
derived from naturally weighted VLA data (after correction for 
primary beam attenuation) by integrating the emission in each 
velocity channel within a 7.'5x 17.'8 aperture centered at c?j2ooo=05 h 
ll m 22.5 s , <5j2ooo = -11°55'01.7". The blue edge of the RX Lep 
profile is blended with Galactic emission. The dotted lines and 
thin solid line show the results of a two-Gaussian fit to the blended 
profile. The spectrum plotted as a dashed histogram was derived 
by summing the emission within irregularly shaped blotches whose 
outer boundaries were defined by 3<r intensity contours. 

gation of the emission and Gaussian line decompositions 
to disentangle the signals. We have verified through mea- 
surements of continuum sources in the field ( £13.3. ID that 
the discrepancy does not result from an uncertainty in 
the amplitude scale of the VLA data. The distinctive 
morphology of the circumstellar H I emission seen in Fig- 
ure 2] also seems to rule out significant contamination 
of the VLA flux density measurements from line-of-sight 
emission. 

4.2.2. Discussion of RX Lep Results 

Based on single-dish mapping and numerical modeling, 
Libert et al. (2008) deduced that RX Lep is surrounded 
by a detached H I shell comprising material from the stel- 
lar wind, augmented by gas swept up from the local ISM. 
They also found indications of elongation of the emission 
along the direction of space motion of the star. Our VLA 
observations of RX Lep confirm these general features of 
the H I distribution. In addition, the wealth of additional 
detail revealed by our imaging data provides important 
new clues on the evolutionary history of the star and its 
interaction with the local environment. 

To first order, the morphological properties of the Hi 
emission surrounding RX Lep appear to be significantly 
influenced by the interaction between the star and its 
local interstellar environment. This is not surprising 
given the moderately high space velocity of the star 
(K;pace ~57 km s _1 ; Table 5). Manifestations of this 
interaction include the elongation of its H I shell, the dis- 
placement of the star from the center of this shell, and 
the extended tail of gas trailing to the south. We now 
discuss each of these properties in turn. 

The north-south elongation of RX Lep's egg-shaped 
shell aligns roughly with the space motion vector that 
we have computed for the star (PAw 25°). This small 
misalignment between the two may be due to a bulk 



ISM flow in the region, directed roughly perpendicular 
to the trajectory of the star. Evidence for similar lo- 
cal flows has recently been inferred from observations of 
other evolved stars (Ueta 2008; Ueta et al. 2010; Wareing 
2012; Menten et al. 2012; see also q4.8.2|) . 

As noted above, RX Lep is displaced by ~ 30" 
(0.022 pc) north of the geometric center of its Hi shell. 
If we adopt half of the FWHM linewidth of the H I spec- 
tral profile as a measure of the outflow velocity (i.e., 
V ou t ~ 2.6 km s _1 ) and a mean shell radius of 0.15 pc, 
the dynamical crossing time for the shell is ^56,000 yr. 
If RX Lep has traveled from the shell center to its 
present position in this amount of time, its implied rate 
of motion in declination would be only ~0.5 mas yr -1 — 
significantly less than its Hipparcos proper motion along 
this direction (56.93 mas yr -1 ). This implies that the 
Hi shell is not "free-floating", but rather that it is be- 
ing dragged along with the moving star. Furthermore, 
this suggests that the displacement of the star from the 
shell center originates primarily from the lower pressure 
encountered in the direction opposite the motion of the 
star compared with in the leading direction. 

As described below, our new data support the pos- 
sibility that RX Lep has only recently begun a new 
phase of mass-loss. If true, this may help to explain an- 
other prominent feature of its shell: the dearth of emis- 
sion along its equatorial region. Villaver et al. (2012) 
point out that for mass-losing stars moving supersoni- 
cally through the ISM, an ensuing pressure drop inte- 
rior to an existing shell that occurs in-between mass- 
loss episodes may trigger instabilities, leading it to break 
apart. Villaver et al. also note that as the wind from 
a new mass loss episode propagates through previously 
ejected material, shock regions are expected to develop. 
We speculate that the break in RX Lep's shell may be 
linked with such shocks or instabilities. 

Extended gaseous wakes resembling cometary-like tails 
have now been detected in Hi emission toward sev- 
eral AGB stars, trailing their motions through space 
(Matthews & Reid 2007; Matthews et al. 2008, 2011b; 
see also H4.8\i . However, none of the tails discovered to 
date has been as highly collimated as the one seen trail- 
ing RX Lep. It is likely that this narrow tail results from 
the confluence of the relatively early evolutionary sta- 
tus of RX Lep compared with other stars imaged in H I 
to date, coupled with its relatively high space velocity 
and the star's location near the Galactic disk (see Ta- 
ble 5). In general, higher stellar space velocities allow 
more efficient deflection of mass loss debris behind the 
star to form a trailing wake (e.g., Villaver et al. 2012), 
while higher local ISM densities (as are expected near the 
Galactic plane) produce greater external pressure, lead- 
ing to greater collimation (e.g., Wareing 2012). We note 
however that in simulations, the adopted cooling function 
and the temperature assumed for the stellar wind will 
also influence the resulting properties of the cometary 
tail. 

A counterpart to the RX Lep tail is also visible in IRAS 
60/im and 100/im images of the star, extending ~ 12' to 
the south (see Figure 9 of Libert et al. 2008). While 
Libert et al. cautioned that this structure could be an 
artifact resulting from the north-south scanning direction 
of the IRAS satellite, our Hi imaging now provides cor- 
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roborating evidence that this feature is likely real. The 
poor angular resolution of the IRAS data preclude a de- 
tailed comparison between the FIR and H I emission, but 
the data are consistent with the two tracers having com- 
parable linear extents. 

To account for their original H I observations of 
RX Lep, Libert et al. (2008) proposed a simple model 
comprising a freely expanding wind with a constant out- 
flow speed and mass loss rate, interior to some radius, r\, 
at which a termination shock occurs owing to the colli- 
sion between the freely expanding stellar wind and the 
surrounding medium. Exterior to T\, the wind sweeps 
up additional matter from its surroundings and is de- 
celerated and compressed as a result of this interaction, 
leading to the formation of a detached shell. 

Incorporating new constraints derived from our VLA 
measurements together with the molecular wind param- 
eters from Libert et al. (2008), we have computed a 
revised version of the Libert et al. model. Details are 
described in Appendix A. We recognize that our model is 
an oversimplification, since it assumes spherical symme- 
try and does not account for the space motion of RX Lep. 
The stellar motion will result in a higher effective pres- 
sure on the leading edge of its H I shell and a pressure in 
the trailing direction that is lower than the average ISM 
pressure. These effects will impart an elongation to the 
H I shell (as well as a reduction of the size of the shell 
in the leading direction) and lead to the formation of 
a trailing wake. Despite these limitations, our model is 
able to provide a satisfactory fit to the globally averaged 
H I line profile of RX Lep (Figure EEh) and supplies a first 
order estimate of the duration of mass loss time history 
for this star. Based on our best-fit model (Table Al), 
we find a mass loss timescale of —90,000 yr. Our derived 
age is consistent with the hypothesis that RX Lep has 
begun its life as a TP-AGB star only relatively recently. 

4.3. Y Ursae Majoris (Y UMa) 

Y Ursae Majoris (Y UMa) is a semi-regular variable 
with an oxygen-rich chemistry. Its spectral type (M7II- 
III) is later than the vast majority of semi-regular vari- 
ables (Dickinson & Dinger 1982). Kiss et al. (1999) 
found this star to have a triply periodic variability with 
a dominant period of 162 days, plus additional modula- 
tions on time scales of 315 and 324 days. (Another star 
in our current sample, the carbon-rich Y CVn, is also a 
triply periodic semi-regular variable; see £)4.5I) . Y UMa 
has an effective temperature of —3100 K based on the 
temperature scale of normal giants of type later than M6 
(Perrin et al. 2004). We adopt a distance to the star of 
385 pc based on the Hipparcos parallax of 2.60±0.59 mas 
(van Leeuwen 2007). 

Y UMa was observed in the CO(2-l) line by Knapp 
et al. (1998). Based on those data, the authors derived 
a stellar systemic velocity of T4.lsr == 19.0±0.5 km s _1 , 
a wind outflow speed of V ou t = 5.4±0.6 km s _1 , and a 
mass-loss rate of M = 2.6 x 1O _7 M yr -1 (scaled to our 
adopted distance). The systemic velocity and outflow 
velocity are in agreement with those derived from the 
thermal SiO u=0 J=2-l and J=3-2 lines detected by 
Gonzalez Delgado et al. (2003). The Hipparcos proper 
motions of —4.95 ± 0.53 mas yr -1 in right ascension and 
2.05±0.50 mas yr -1 in declination imply a peculiar space 



motion of 19.2 km s 1 along a position angle of 57°1 
(Table 5). 

4.3.1. VLA Results for Y UMa 

We present H I channel maps from our naturally 
weighted VLA data cube in Figure [6] Emission is de- 
tected at the stellar position in the channel centered at 
Ulsr = 19-0 km s" 1 (corresponding to the stellar sys- 
temic velocity) and in several surrounding channels. In- 
deed, we detect Hi emission at > 3cr that appears to be 
associated with the CSE over the velocity range 13.8 to 
22.9 km s _1 . The velocity extent of the Hi emission is 
thus nearly identical to that of the CO(2-l) emission de- 
tected by Knapp et al. (1998). Redward of the stellar 
velocity (Vlsr > 19-6 km s -1 ), the detected emission is 
weak and unresolved and lies at or near the stellar po- 
sition. However, for velocities Vlsr < 19-0 km s _1 , the 
emission is clearly extended, and there are hints of a bi- 
furcation in the emission structure. This is most clearly 
seen in the 15.8 km s _1 and 16.4 km s^ 1 velocity chan- 
nels. 

In Figure [7] we present an H I total intensity map of 

Y UMa. This map was produced from the robust +1 
weighted data (Table 4) to highlight the clumpy nature of 
the emission and the existence of a "notch" in the CSE to 
the northeast of the star, along a position angle of ~ 60°. 
This position angle closely matches that of the direction 
of space motion of the star (Table 5). An analogous 
feature, also along the direction of space motion, is seen 
in the CSE of Y CVn (Figure US] gXp) . One can also 
see from Figure [7| that the peak H I column density in 
the CSE of Y UMa is displaced by - 20" (-0.037 pc) to 
the northeast of the stellar position. 

The morphology of the H I emission surrounding 

Y UMa suggests that the emission in the northwestern 
part of the CSE has been severely distorted by its in- 
teraction with the surrounding medium or other exter- 
nal forces. The result is a CSE with a distinctly non- 
spherical shape, including a small plume of gas extend- 
ing from the western edge of the H I nebula. Despite this 
distortion, we measure nearly equal H I masses on either 
side of the bisector running along the position angle of 
the space motion. Including the extended plume, the 
maximum Hi extent of Y UMa is ~6f8 (0.76 pc). 

Figure [8] shows an H I velocity field for Y UMa, again 
derived from the robust +1 data. The Hi velocity field 
is clearly dominated by emission blueward of the stel- 
lar systemic velocity. Interestingly, the velocity pattern 
exhibits a 180-degree rotational symmetry, but with the 
symmetry point displaced by — V northeast of the stel- 
lar position. This displacement is along the space motion 
vector of the star. One can also see in Figure |8] evidence 
for systematic velocity gradients across the CSE; along 
PA« 70° and 250°, the emission becomes systematically 
redshifted with increasing distance from the star, while 
along PAw 110° and « 350°, the emission becomes sys- 
tematically blueshifted with increasing distance. These 
angles match the bifurcation visible in the channel im- 
ages in Figure El Y UMa is the only star imaged in H I to 
date to show these types of features in the velocity field, 
and their origin and interpretation are presently unclear. 

A spatially integrated H I spectrum of Y UMa is pre- 
sented in Figure [9] As expected from the emission dis- 
tribution seen in Figure [71 the peak integrated flux den- 
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Fig. 6. — Hi channel maps bracketing the systemic velocity of Y UMa, derived from naturally weighted VLA data. The spatial resolution 
s ~ 72" X 58". A star symbol indicates the stellar position. Contour levels are (-4.2,-3,3,4.2,6,8. 4)xl. 7 mjy beam -1 . The lowest contour 
s ~ 3a. The field-of-view shown is comparable to the VLA primary beam. 



sity occurs blueward of the stellar systemic velocity. The 
spectrum has a two-component structure similar to what 
was reported previously by Gerard & Le Bertre (2006) 
based on NRT spectroscopy. Using a two-component 
Gaussian fit to the VLA spectrum, we find a narrow 
component with a peak flux density of 0.10±0.01 Jy 
and a FWHM line width of 3.2±0.4 kms" 1 , centered 
at Vlsr = 16-4 ± 0.2 km s" 1 , and a weaker, broad 
component centered at 17. 2 ±0.8 km s^ 1 with a peak 
flux density 0.03±0.01 Jy and a FWHM line width of 
9.2±2.3 km s _1 . Both components are blueshiftcd rel- 
ative to the systematic velocity of the star as deter- 
mined from CO observations, although the broad com- 
ponent matches the systemic velocity determined from 
near infrared lines by Lebzelter & Hinkle (2002) (14 = 
17.7±0.5 km s" 1 ). 

The peak Hi flux density measured with the VLA is 
~40% higher than the original NRT measurement by 
Gerard & Le Bertre (2006). However, a new analysis us- 
ing higher quality NRT data now shows good agreement 
between the VLA and NRT measurements, although the 
NRT data hint that the source size may be more ex- 



tended than we observe with the VLA (up to 16'; Gerard 
& Le Bertre, in prep.). 

Based on our Gaussian fit to the spectrum in Figure [9j 
we measure an integrated H I flux density associated with 
the CSE of Y UMa of / S ffl dv = 0.63±0.13 Jy km s _1 , 
translating to a total Hi mass A/hi ~ 0.022 M . Al- 
though the double Gaussian provides a good fit to the 
line profile, the peak of the composite Gaussian profile is 
shifted by ^0.6 km s _1 to the blue of the observed line 
peak, and there are hints of residual emission in the blue 
wing of the line. 

4.3.2. Discussion of Y UMa Results 

From Figure |5] one sees that the emission giving rise 
to the broader of the two Hi line components in Fig- 
ure [9] arises from spatially compact emission lying at or 
near the position of the star, while the narrower veloc- 
ity component comes from material that is more spa- 
tially extended. This suggests that the broad line com- 
ponent may be linked with the freely expanding wind 
detected at small radii through its CO and SiO emis- 
sion. This hypothesis is consistent with the respective 
widths of the Hi and molecular lines. If we approxi- 
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Fig. 7. — Hi total intensity map of Y UMa derived from data with robust +1 weighting. The spatial resolution is ~ 61" X 52". The 
map was constructed from emission spanning LSR velocities 13.8 to 22.9 km s — 1 . Contour levels are (l,1.4,2,2.8,4,5.6,8)x4.8 Jy beam -1 
m s -1 and the intensity levels are to 35 Jy beam -1 m s — . To minimize noise in the map, data at a given point were blanked if they did 
not exceed a 1.8<r threshold after smoothing by a factor of three spatially and spectrally. A star symbol indicates the stellar position and 
white line indicates the direction of space motion. The field-of-view shown is comparable to the VLA primary beam. 



mate the molecular linewidth as twice the CO outflow 
velocity (~10.8±0.6 km s _1 ), this matches the FWHM 
of the broad Hi line, although the line centroids differ 
by ^2-3 km s -1 . Higher resolution Hi imaging together 
with CO imaging of Y UMa would be of interest for more 
thoroughly investigating this link. 

Y UMa was reported by Young et al. (1993a) to be an 
extended infrared source with a diameter of T.Q at 60/im. 
The extended nature of the FIR emission is confirmed by 
an archival ISO 65/im image shown in Figure [TU1 To our 
knowledge, these ISO data have not been published pre- 
viously. The observations were executed on 1996 October 
18 as part of program 33700231 (Principal Investigator: 
J. Gurtler). We have downloaded the oversampled map 
as a "Standard Processed Data" FITS file from the ISO 
data archived The original image had a field-of-view of 
687" x 687" and a pixel size of 15"0. For display pur- 

8 The archive is h osted by the European Space Agency at 
http:/ /iso.csac.csa.int/ida 



poses, the image was resampled to a pixel size of 10" to 
match the VLA data. 

While the FIR emission appears to have a roughly com- 
parable spatial extent as the Hi emission detected with 
the VLA, the two do not correspond closely in detail. 
For example, the FIR emission appears elongated along 
a northeast /southwest direction, extending further to the 
southwest than the H I emission. Conversely, we see that 
there is little or no FIR emission underlying the north- 
western portion of the H I nebula, including the ridge 
of enhanced Hi column density and the weaker plume 
extending from the northern edge. Since in principle, 
whatever process is responsible for the distortion on the 
northwestern side of the CSE and creation of the "plume" 
observed in H I should have acted on both the H I and the 
FIR-emitting (presumably dusty) material, this suggests 
that either the atomic hydrogen is intrinsically more ex- 
tended than the dust or that dust grains were destroyed 
during this process. The absence of any strong correla- 




Fig. 8. — Hi velocity field of Y UMa derived from naturally weighted data. The spatial resolution is ~ 72" X 58". The intensity scale 
indicates LSR radial velocity in km s _1 . 
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Fig. 9. — Spatially integrated Hi spectrum of Y UMa. The 
spectrum was derived from naturally weighted VLA data (after 
correction for primary bean attenuation) by integrating the emis- 
sion in each velocity channel within a 8.'5 X 7.'0 aperture centered at 
aj2000=12 h 40 m 20.1 s , <5j2000 =55° 51' 27.6". The dotted curves 
and thin solid curve show, respectively, a two-component Gaussian 
fit to the line profile and the sum of the two components. The 
vertical bar indicates the stellar systemic velocity as determined 
from CO observations. 

tions between the atomic gas and the FIR emission may 
also imply that the dust grains are not uniformly mixed 
throughout the CSE. Simulations by van Marie et al. 
(2011) have shown that such non- uniform mixing is ex- 
pected for the largest CSE grains (with radii ^0.105/im). 



As in the case of RX Lep ( HA.2. 2[) . we have derived a 
simple model for the H I emission of Y UMa under the 
assumption that it arises from an outflowing wind that 
has been abruptly slowed at a termination shock where 
it meets the surrounding ISM, resulting in the formation 
of a detached shell (see Appendix A and Table Al for 
details). Despite the observed distortions of the Y UMa 
H I shell, our simple model provides a reasonable fit to the 
global (spatially averaged) Hi line profile, including its 
two-component structure (Figure [2"8"b ). Using this model 
we estimate a total duration of the mass loss history of 
the star of ~ 1.2 x 10 5 yr. 

4.4. R Pegasi (R Peg) 

R Pegasi (R Peg) is an oxygen-rich Mira variable with 
a mean spectral type M7e and a pulsation period of 378 
days (Jura & Kleinmann 1992; Whitelock et al. 2008). 
The effective temperature of R Peg is known to vary 
significantly throughout its pulsation cycle; van Belle et 
al. (1996) measured T e g = 2333±100 K at minimum and 
2881±153 K at maximum. 

We adopt a distance to R Peg of 400 pc based 
on the period-luminosity relation for Miras derived by 
Jura & Kleinmann (1992), although we note that the 
study of Whitelock et al. (2008) yielded a value -25% 
higher (490 pc). The Hipparcos parallax for this star 
(2.61±2.03 mas yr -1 ; Knapp et al. 2003) is too uncer- 
tain to provide any additional meaningful constraints on 
the distance 

Based on the CO(2-l) observations of Winters et 
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Fig. 10. — Hi total intensity contours (same as in FigureO overlaid on an ISO 65/xm image of Y UMa. The brightness units of the ISO 
image are MJy steradian - 1 . 



al. (2003), we assume a stellar systemic veloc- 
ity of 14 ! lsr=24.0 km s _1 , an outflow velocity of 

V ou t=5.5 km s _1 , and a mass loss rate of M ~ 5.3 x 
10~ 7 M yr -1 (scaled to our adopted distance). Us- 
ing the Hipparcos proper motions from van Leeuwen 
(2007) of 12.63±1.74 mas yr _1 in right ascension and 
—9.71 ±1.67 mas yr _1 in declination, we derive a peculiar 
space motion of T4 pacc =26.3 km s~ 2 along a PA=147°6 
(Table 5). 

R Peg was previously detected in Hi by Gerard 
& Le Bertre (2006), who found evidence for a two- 
component line profile comprising a strong, narrow, 
roughly Gaussian component superposed on a broader 
pedestal. They found a velocity offset of ^2 km s _1 be- 
tween the two components, with the two bracketing the 
systemic velocity of the star as determined from CO ob- 
servations. 

4.4.1. VLA Results for R Peg 

Hi channel maps surrounding the stellar systemic ve- 
locity of R Peg are presented in Figure [TTJ The mor- 
phology of the H I emission seen in these channel images 
is unlike any other star that has been imaged in H I to 
date. H I emission is clearly detected near the stellar po- 
sition in several velocity channels bracketing V*. How- 
ever, at VlSR=25.1 km s _1 , the emission peak is offset 
by roughly an arcminute (one beam width) to the south- 
east of the stellar position. In the three velocity chan- 
nels from 23.2 to 24.5 km s _1 , multiple emission peaks 
are present near the stellar position, but again the peak 



emission in each channel is offset to the southwest of the 
star, with a relative dearth of emission at the stellar po- 
sition itself. The morphology in these channels is some- 
what reminiscent of a bipolar flow (but see below). This 
bifurcation is even more evident at Vlsr=22.5 km s _1 , 
but here the far reaches of the two lobes curve to form 
a horseshoe shape, with an extent of > 15' (1.7 pc). 
This "horseshoe" is also delineated by multiple emission 
clumps at VLsr=21.9 km s" 1 . Finally, from Vlsr=20.0 
to 21.3 km s _1 , the emission clumps along the horseshoe 
decrease in number and peak column density while their 
projected distances from the star increase. 

We present an Hi total intensity map for R Peg in 
Figure This map confirms that the peak H I column 
density toward R Peg does not coincide with the stellar 
position, but instead lies at a projected distance of ~ V 
(0.1 pc) southeast of the star. 

A global H I spectrum of R Peg is shown in Figure [T3l 
There are hints of a slight asymmetry in the line pro- 
file, as well as a line peak that is broader and flatter 
than the other stars in the current sample. Nonetheless, 
a single Gaussian (dotted line on Figure [T3| provides a 
reasonable fit to the data for the purpose of characteriz- 
ing the global line parameters. Based on such a fit, we 
measure a peak Hi flux density of 0.16±0.01 Jy, in ex- 
cellent agreement with Gerard & Le Bertre (2006). We 
also find a FWHM line width of 4.35±0.39 km s~\ cen- 
tered at VLsr=22.4±0.2 km s" 1 . The Hi line centroid is 
blueshifted by ~1 km s _1 relative to the stellar systemic 
velocity determined from CO observations. From our 
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Fig. 11. — H I channel maps bracketing the systemic velocity of R Peg, derived from naturally weighted VLA data. The spatial resolution 
is ~ 63" X 57". A star symbol indicates the stellar position. Contour levels are (-4.2[absent],-3,3,4.2,6,8.4)xl.7 mjy beam -1 . The lowest 
contour is ~ 3o\ The field-of-view shown is comparable to the VLA primary beam. 



Gaussian fit we also calculate an integrated H I flux den- 
sity of J Sm dv =0.76±0.09 Jy km s -1 , translating to an 
H I mass of Mm ~ 0.029 M Q . After scaling to account for 
our different adopted distances, the latter agrees closely 
with the value of 0.03 M Q previously reported by Gerard 
& Le Bertre (2006). 

Gerard & Le Bertre (2006) also reported a second 
weaker, broader line component centered on the stel- 
lar systemic velocity and with a velocity extent com- 
parable to that of the CO(2-l) line emission (FWHM 
~10 km s _1 ). New NRT measurements (Gerard & 
Le Bertre, in prep.) reaffirm the presence of this broad 
component with a peak flux density of 0.006±0.003 Jy, 
too weak to have been detected in the current VLA data. 



4.4.2. Discussion of R Peg Results 

Above we noted the presence of an offset of ~0.1 pc 
to the southeast between the peak Hi column density 
observed near R Peg and the stellar position. The di- 
rection of this offset is within ~ 10° of position angle of 
the space velocity vector that we have computed for the 
star (PA« 148°; Table 5), suggesting that one possible 
explanation is that there is an accumulation of mate- 
rial in front of the star caused by the "snowplow effect" 
(Isaacman 1979). In this picture, the gas to the north- 
west of the star could comprise a trailing wake of ma- 
terial, ram pressure-stripped through the interaction be- 
tween the CSE and the local ISM (see £|4.2.2j) . Consistent 
with this hypothesis, over the velocity channels spanning 
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Fig. 12. — Hi total intensity map of R Peg derived from naturally weighted data. A star symbol indicates the stellar position and the red 
line indicates the direction of space motion. The spatial resolution is ~ 63" X 57". The map was constructed from emission spanning LSR 
velocities 20.0 to 25.1 km s -1 . Contour levels are (1,1.4,2,2.8,4,5.6,8) x4 Jy beam -1 m s and the intensity range is to 43 Jy beam -1 
m s — . To minimize noise in the map, data at a given point were blanked if they did not exceed a 2a threshold after smoothing by a factor 
of three spatially and spectrally. 



25.1 km s -1 to 22.5 km s , the emission in the putative 
wake shows a systematic increase in projected distance 
from R Peg with decreasing LSR velocity, as would be 
expected for debris that is being slowly decelerated by 
dynamical friction from the local ISM (see Matthews et 
ai. 2008, 2011b; Raga & Canto 2008). A related possi- 
bility is that the lobes of H I emission near the star are 
the remnants of a wake of debris shed during a previous 
cycle of mass-loss, and that the material has since drifted 
relative to the stellar position. This could also help to 
explain the small blueshift of the H I line centroid relative 
to the systemic velocity determined from molecular line 
observations visible in Figure H3l However, an ISM flow 
originating from a northerly direction would seem to be 
required to account for the direction of the offset between 
the star and the H I emission, as well as the curving sides 
of the horseshoe, which are oriented roughly perpendic- 
ular to the direction of the stellar space motion. 

Despite the bifurcated appearance of the emission in 
several of the Hi channel maps (Figure ITT]) , we con- 



sider the possibility that the Hi surrounding R Peg re- 
sults from some type of bipolar outflow to be unlikely. 
While bipolar flows arc known in the case of a number of 
AGB stars, with the exception of the high- velocity flows 
(~100 km s -1 ) arising from certain known or suspected 
binaries (Imai et al. 2002; Sahai et al. 2003), these 
outflows do not tend to be strongly collimated (e.g., In- 
omata et al. 2007; Libert et al. 2010; Castro-Carrizo 
et al. 2010). Further, in the case of R Peg, the two 
"lobes" would have to be oriented close to the plane of 
the sky to account for their similar radial velocities. Fi- 
nally, the asymmetry of the velocity-integrated emission 
relative to the stellar position and the elongation of the 
emission parallel to the space motion vector (Figure H2|) 
are both effects expected as a result of interaction be- 
tween the ejecta of a moving star and the surrounding 
ISM. 

4.5. Y Canum Venaticorum (Y CVn) 
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Fig. 13. — Spatially integrated Hi spectrum of R Peg. The spec- 
trum was derived from naturally weighted VLA data (after cor- 
rection for primary beam attenuation) by integrating the emission 
in each velocity channel within a 15.'5x8.'3 aperture centered at 
aj2000=23 h 06 m 18.49 s , 5 J2 ooo = 10°31'06.0". The dotted line 
shows a Gaussian fit to the line profile. The vertical bar indicates 
the stellar systemic velocity as determined from CO measurements. 

Y Canum Venaticorum (Y CVn) is a rare example of a 
J-type carbon star. Its exact evolutionary stage is poorly 
known (e.g., Dominy 1984; Lambert et al. 1986; Lorenz- 
Martins 1996). The absence of technetium absorption 
(Little et al. 1987) and the lack of enhancement in s- 
process elements (Utsumi 1985) both suggest that the 
star has not yet undergone a thermal pulse. Some au- 
thors have even suggested that Y CVn is still on the 
red giant branch and that its carbon-rich composition 
may have resulted from a core He flash (Dominy 1984), 
although the high derived mass loss rate from the star 
appears to preclude this possibility and place the star on 
the AGB (e.g., Izumiura et al. 1996). The bolometric 
luminosity of Y CVn also places it on the AGB; Libert 
et al. (2007) derived a value of 9652L (scaled to our 
adopted distance; see below) using the bolometric cor- 
rection of Le Bertre et al. (2001). 

Y CVn belongs to the class of semi-regular variables, 
having a dominant pulsation period of ^160 days and 
secondary periods of 273 and 3000 days (Kiss et al. 
1999). (Another star in our sample, the oxygen-rich 
Y UMa, is also a triply periodic semi-regular variable; 
see i)4.3p . Lambert et al. (1986) derived an effective 
temperature for Y CVn of ^2730 K. 

From CO(3-2) observations, Knapp et al. (1998) de- 
rived a stellar systemic velocity for Y CVn of T4,lsr 
21.1 ■ ().!) km s -1 , an outflow velocity for its wind of 
Vo U t=7.8±1.3 km s _1 , and a mass loss rate M = 1.7 x 
10~ 7 Mq yr" 1 (scaled to our adopted distance). We note 
however that the observations of Jura et al. (1988) im- 
ply a large 13 C abundance in the star, thus the mass-loss 
rate derived from 12 CO lines may be underestimated by 
a factor of ~1.4 (see Libert et al. 2007). 

The distance to Y CVn is uncertain, with values in 
the literature differing by more than ±60 pc. Here we 
adopt a value of 272 pc based on the re-analysis of the 
Hipparcos parallax measurements by Knapp et al. (2003; 
3.68±0.83 mas). The proper motion is —1.40 ±0.31 mas 
yr _1 in right ascension and 13.24±0.26 mas yr _1 in decli- 



nation (van Leeuwen 2007), which translates to a peculiar 
space velocity of 14 pace =31.0 km s _1 at PA=36°8. 

Y CVn was previously detected in H I with the NRT by 
Lc Ber tre & Gerard (2004) and Libert et al. (2007) (see 
£ ]4.5.2l for a discussion). Y CVn has also long been known 
to be surrounded by extended FIR emission (Young et 
al. 1993a; Izumiura et al. 1996; Geise 2011; Cox et al. 
2012a). The ISO 90/xm data of Izumiura et al. show 
a central infrared source surrounded by a quasi-circular 
shell with an inner radius r; n ^2'.8 and an outer radius 
»"out ~5^1. This shell is displaced slightly relative to the 
central infrared source. Izumiura et al. interpreted the 
properties of the FIR emission as indicating that Y CVn 
has undergone a decrease in mass loss rate of two orders 
of magnitude compared to the value at the time of the 
dust shell formation. However, Libert et al. proposed 
an alternative interpretation, namely that the shell is a 
several hundred thousand year-old structure formed by 
the slowing down of the stellar wind by the surrounding 
matter. This does not rule out that changes in the mass 
loss rate over time have occurred in Y CVn, but it re- 
moves the requirement of an abrupt drop in M posterior 
to the formation of shell in order to explain the central 
hole seen in the FIR observations. As described below, 
our new data strengthen the interpretation put forth by 
Libert et al. and permit a more detailed comparison with 
the ISO FIR data. 

4.5.1. VLA Results for Y CVn 

Hi channel maps for Y CVn derived from a tapered, 
naturally weighted version of our VLA data cube are 
presented in Figure [14] H I emission associated with the 
CSE of Y CVn is detected at LSR velocities ranging from 
15.3 to 25.6 km s _1 . We show a tapered version of our 
data here despite their lower spatial resolution, as they 
better highlight the faint emission features present near 
the stellar position in the velocity channels from 15.3 to 
17.9 km s _1 and from 24.3 to 25.6 km s _1 . At the ve- 
locities closest to the stellar systemic velocity, the emis- 
sion appears extended and clumpy and distributed over 
a shell-like structure, while a relative dearth of emission 
is notable directly toward the stellar position. 

An H I total intensity map of Y CVn derived from natu- 
rally weighted data is presented in Figure [15] The stellar 
position is offset ~ 45" to the northeast relative to the 
geometric center of the H I emission distribution. The 
H I nebula is roughly symmetric about the position angle 
of this offset (~42°), which corresponds to the direction 
of space motion of the star (Table 5). In contrast to the 
quasi-spherical appearance of the CSE in previous FIR 
images (Izumiura et al. 1996), the integrated Hi distri- 
bution exhibits a more boxy shape, with the extent of the 
northeastern edge of the nebula being somewhat broader 
(~ 9') than the southwest side (~ 7'). This boxiness is 
also visible in the Spitzer 70/im image of Geise (2011). 

The H I distribution of Y CVn is quite clumpy down to 
the resolution limit of our data, with two of the brightest 
clumps lying at the southern and western "corners" of 
the nebula. Several additional compact, quasi-spherical 
clumps of gas are also detected outside the main emission 
distribution and appear to be connected with the CSE 
by faint (~ 3er) tendrils of gas. The nature of these 
structures is unclear, but they appear to be circumstellar 
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Fig. 14. — Hi channel maps bracketing the LSR velocity of Y CVn, derived from tapered, naturally weighted VLA data. The velocities 
shown correspond to the values over which CO(3-2) emission was detected by Knapp et al. 1998. A star symbol indicates the stellar 
position. The spatial resolution is ~ 107" X 83". Contour levels arc (-3.5,-2.5,2.5,3.5,5,7.1,10.0,14.1) x2. 6 mjy beam -1 . The lowest 
contour is ~ 2.5<r. The greyscale intensity levels are 2.7 to 30.0 mjy beam -1 . 
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Fig. 15. — H I total intensity map of Y CVn derived from naturally weighted VLA data. The spatial resolution is ~ 61" X 49". The map 
was constructed from emission spanning LSR velocities 17.8 to 23.0 km s _1 . Contour levels are (0.5,0.7,1)1-4,2,2.8,4,5.6,8) X 13 Jy beam 
m s~ 1 . To minimize noise in the map, data at a given point were blanked if they did not exceed a 2<r threshold after spatially smoothing 
the data by a factor of three. Intensity levels are 4 to 68 Jy m s _1 . A star symbol indicates the stellar position, and the red line indicates 
the direction of space motion. 
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Fig. 16. — Spatially integrated Hi spectra of Y CVn, derived 
from naturally weighted VLA data (after correction for primary 
beam attenuation). Upper panel: Spectrum obtained by integrat- 
ing the emission in each velocity channel within a 12.'0xl2.'2 aper- 
ture centered at Q J20 oo = 12 h 45 m 01.2 s , <5 J20 oo = 45°25'39.9". The 
dotted line shows a Gaussian fit to the spectrum and the vertical 
bar indicates the stellar systemic velocity as determined from CO 
observations. Lower panel: Spectrum integrated within a circular 
aperture of radius 2.'8 centered on the star. The dotted lines show 
a two-component Gaussian fit and the thick solid line shows their 
sum (see text for details). 

in nature. 

The depression in the total H I column density toward 
the stellar position seen in the individual H I channel 
maps is also apparent in the H I total intensity map. In 
addition, the Hi total intensity map reveals a "notch" 
extending from the star to the northeastern edge of the 
nebula along PAw 42° where relatively little H I emis- 
sion is detected. The nature of this feature is discussed 
further below. 

The spatially integrated H I spectrum of Y CVn (Fig- 
ure [16]) displays a narrow, roughly Gaussian shape, simi- 
lar to what we observe for V1942 Sgr ( 34.6.2[) . Based on 
a single Gaussian fit (dotted line on Figure [TBI we mea- 
sure a peak flux density of 0.79±0.02 Jy and a FWHM 
linewidth of 3.2±0.1 km s _1 centered at an LSR veloc- 
ity of 20.65±0.05 km s _1 . From this fit we derive an 
integrated Hi flux density of / Shi dv =2.68±0.11 Jy 
km s _1 , translating to an Hi mass A/hi ~ 0.047 M Q . 

The H I mass and global line width that we derive from 



the VLA data are in excellent agreement with those de- 
rived previously from NRT observations by Le Bertre & 
Gerard (2004), although our integrated Hi line flux is 
~ 25% smaller than reported by Libert et al. (2007). 
Further comparisons with prior NRT results are de- 
scribed below. 

4.5.2. Discussion of Y CVn Results 

Above we drew attention to the presence of a "notch" 
(depression in the H I column density) in the shell of 
Y CVn, lying ahead of the star and running parallel to 
its direction of space motion. A similar feature was also 
seen in the CSE of Y UMa ( 34.3.20 . In both cases, a 
corresponding depression is also present in the FIR light 
(Figure [10] and Figure |T7l discussed below) . The origin 
of these notches is unclear, but one possibility is that 
they are linked to instabilities along the leading edge of 
the CSE. For example, in the hydrodynamic simulations 
of mass-losing stars moving through the ISM performed 
by Wareing et al. (2007a, b) and Villaver et al. (2012), 
instabilities may develop under certain conditions that 
lead to a bending of the leading edge of the CSE (bow 
shock) back toward the star. Unfortunately, none of the 
currently published simulations are closely matched to 
the current age, mass-loss parameters, and space motion 
parameters of Y CVn and Y UMa. 

Another interesting phenomenon seen in the Hi im- 
ages of Y CVn is the presence of multiple compact knots 
of emission lying just exterior to the CSE at several lo- 
cations along both its leading and trailing edges. These 
knots are detected at high significance and appear to be 
linked to the Y CVn shell by tendrils of gas, suggesting 
that they are quite likely to be circumstellar in nature. 
One possibility is that these knots may arise from ther- 
mal instabilities (e.g., Zucker & Soker 1993) along the 
interface between the Y CVn shell and the surrounding 
ISM. 

In Figure [T7] we present Hi contours derived from the 
robust +1 version of our VLA data overlaid on the ISO 
90/im data from Izumiura et al. (1996). We downloaded 
the archival data for Y CVn as described in 34.3.21 The 
original ISO image has a pixel size of 15"0 but was re- 
binned to 10" pixels for display purposes. 

The shape and extent of the Hi and FIR emitting 
shells are quite similar along the northwest and south- 
west edges. However, along the southeast portion of the 
CSE, the Hi contours extend outside the FIR emission 
detected by ISO. At the same time, an absence of Hi 
emission is apparent along the northeast (leading) side 
of the FIR-emitting shell. The gas along this edge may 
have been ionized by shocks owing to the interaction with 
the ISM. However, there is no obvious enhancement in 
FIR surface brightness or change in morphology along 
the leading edge of the shell that would indicate the pres- 
ence of a strong bow shock capable of ionizing hydrogen. 
Similarly, no clear evidence of a FIR bow shock is seen 
in the more recent, higher resolution observations from 
Herschel (Cox et al. 2012a). We have also examined 
archival Galaxy Evolution Explorer (GALEX) images of 
the region around the star to search for FUV emission 
from a possible bow shock as has been detected in the 
case of at least three evolved stars (Martin et al. 2007; 
Sahai & Chronopoulos 2010; Le Bertre et al. 2012), but 
no extended FUV emission is visible. 
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Fig. 17. — Hi total intensity contours of Y CVn derived from VLA data with robust +1 weighting, overlaid on an ISO 90/xm image of the 
star from Izumiura et al. 1996. The Hi data have a resolution of ~ 54" X 44". The contour levels are (1,1.4,2,2.8,4,5.6,8) X 13 Jy beam -1 
m s" 1 . The ISO image has brightness units of MJy sr — 1 . 



Izumiura et al. (1996) noted that the western part 
of the FIR shell of Y CVn has a slightly lower surface 
brightness in both 90/im and 160/xm images compared 
with the eastern side. Similarly, we find that the Hi 
emission is slightly enhanced on the eastern side of the 
shell compared with the western side (cf. Figures [15] & 
[T7|) . After applying a multiplicative correction of 1.4 to 
account for the mass of helium, the mass we derive for 
the circumstellar shell of Y CVn (w 0.066 Mq) is consis- 
tent with the value derived previously from FIR observa- 
tions (Af s hoii ~0.047 to 0.16 Mq, scaled to our adopted 
distance; Izumiura et al. 1996). Since the latter determi- 
nation depends on the adopted dust-to-gas ratio, implicit 
to this agreement is that the dust-to-gas ratio of the shell 
is typical of other carbon-rich stars (~ 4.5 x 10 -3 ; Jura 
1986). It also assumes that the gas in the shell is pre- 
dominantly atomic, although as noted by Libert et al. 
(2007), the temperature and density conditions within 
Y CVn's shell may be conducive to molecule formation 
through a non-equilibrium chemistry. 

Our VLA imaging confirms that the narrow H I line 
component detected by Le Bertre & Gerard (2004) and 



Libert et al. (2007) is closely linked with material dis- 
tributed within the extended FIR-emitting shell sur- 
rounding Y CVn. However, both of these studies also 
reported a second broad (FWHM^16 km s^ 1 ) line com- 
ponent that they attributed to a counterpart to the freely 
expanding wind traced by molecular lines at small radii. 
This broad component is not evident in the global Hi 
profile presented in Figure 1161 However, we tentatively 
identify a counterpart to this component in Figure [16] 
(lower panel), where we plot a spectrum derived from 
only the emission confined within a radius of 2'.8. (This 
radius was adopted based on a model of the shell pre- 
sented in Appendix A; see also below). We find that a 
single Gaussian fit to the spectrum in the lower panel 
of Figure [16] leaves statistically significant residuals in 
the line wings. Adopting instead a two-Gaussian fit, we 
find a narrow component with -Fpcak = 0.35 ± 0.02 Jy 
and a FWHM line width of 3.0 ±0.1 km s" 1 , centered at 
Vlsr = 20.68±0.04 km s -1 (i.e., consistent in position 
and width with the fit to the global spectrum in the up- 
per panel of Figure 1161) and a broad component centered 
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at 20.6±0.5 km s" 1 with F pcak = 0.04 ± 0.02 Jy and a 
FWHM line width of 10.0±1.8 km s _1 . The line width 
we derive for the broader component is smaller than the 
values reported by Le Bertre & Gerard (14.3 km s _1 ) and 
by Libert et al. (~16 km s _1 ). However, Libert et al. 
fitted a rectangular rather than a Gaussian line profile, 
and moreover, the single-dish determinations are sensi- 
tive to uncertainties in baseline fitting. Thus account- 
ing for systematic uncertainties, it appears that all three 
measurements of the broad Hi component are roughly 
consistent with the CO (3-2) line, for which Knapp et 
al. f 998 measured a FWHM~12 km s _1 and derived 
V ou t = 7.8±1.3 km s _1 . This supports the suggestion 
that this material is linked with Y CVn's freely expand- 
ing wind. Corroborating evidence is also seen in our Hi 
channel maps (Figure fT4|) . where weak emission (> 2.5a) 
is detected near the stellar position at velocities offset 
by ~ ±4 to 6 km s _1 from the stellar systemic veloc- 
ity. However, we note that there is evidence that the H I 
emission at these velocity extremes is spatially extended 
on scales > V and thus must be extended well beyond 
the outer boundary of the detected CO(f-O) emission 
(r ~6"5) and the molecular dissociation radius (r ~ 20"; 
Neri et al. 1998). Additionally, the peaks of this emis- 
sion near the velocity extrema in some channels appear 
offset from the stellar position by ~ 1'. 

Our new VLA measurements allow us to provide im- 
proved constraints on the numerical model of the H I shell 
of Y CVn presented by Libert et al. (2007) (see Ap- 
pendix A for details) . In Figure l2"5b we show our model 
spectrum overplotted on the global VLA H I spectrum of 
Y CVn. Based on this model, we estimate an age for the 
H I shell of ~ 4.5 x 10 5 yr. This value is consistent with 
the previous determination of Libert et al., but we note 
that we have adopted a larger distance to the star than 
those authors as well as slightly different model parame- 
ters. Our value is more than 6 times larger than the age 
for the circumstellar shell previously derived by Young 
et al. (1993b) based on FIR data. 

4.6. V1942 Sagittarii (V1942 Sgr) 

V1942 Sagittarii (V1942 Sgr) has been previously clas- 
sified an N-type carbon star and a long period irregular 
(Lb) variable. However, Miller et al. (2012) recently re- 
ported evidence that V1942 Sgr belongs to the class of 
DY Persei-type stars. DY Persei stars, only a handful 
of which have been identified to date, are thought to be 
either a cooler subclass of R Coronae Borealis stars (see 
Clayton 2012) or normal carbon stars that have experi- 
enced ejection events (Tisserand et al. 2009). DY Persei- 
type stars generally exhibit C-rich spectra similar to re- 
type carbon stars and show detectable 13 C. They are 
also characterized by large-amplitude, irregular varia- 
tions in brightness (>1.5 magnitudes). In the case of 
V1942 Sgr, Miller et al. reported evidence for multiple 
periodicities in the light curve variations (120, 175, and 
221 days) that appear to change on time scales of 1-2 
years. However, Soszyhski et al. (2009) found no evi- 
dence for a marked discontinuity between the variations 
observed among DY Persei stars and other carbon-rich 
long-period variables. 

V1942 Sgr has an effective temperature T e g w 2960 K 
(Bergeat et al. 2001) and a bolometric luminosity of 



52OOL0 (Libert et al. 2010a). Its properties are believed 
to place V1942 Sgr on the TP-AGB (e.g., Bergeat et al. 
2002), although its new classification as a DY Persei- 
type star suggests that its evolutionary status and/or 
history could be somewhat different from the other car- 
bon stars previously studied in H 1. We adopt a distance 
to V1942 Sgr of 535 pc based on the Hipparcos parallax 
of 1.87±0.51 mas (van Leeuwen 2007). 

IRAS 60/im data show V1942 Sgr to be surrounded by 
an extended infrared shell of radius 3'.2 (~0.5 pc; Young 
et al. 1993a). V1942 Sgr exhibits a featureless mid- 
infrared (8-22/im) spectrum, consistent with the absence 
of current mass loss (Olnon & Raimond 1986). However, 
the changes in its variability and its DY Persei classifi- 
cation suggest that regular ejection events may be tak- 
ing place (Miller et al. 2012). Furthermore, Libert al. 
(2010a) found evidence of ongoing mass loss from CO 
observations. 

The CO(1-0) and CO(2-l) profiles of V1942 Sgr ob- 
tained by Libert et al. (2010a) exhibit a two-component 
structure comprising broad and narrow components. In 
several cases, such line profiles have been found to be 
the hallmark of bipolar outflows in AGB stars (e.g., Ka- 
hane & Jura 1996; Bergman et al. 2000; Josselin et 
al. 2000; Libert et al. 2010b; C astro- Carrizo et al. 
2010), although CO imaging observations are needed to 
confirm this interpretation for the case of V1942 Sgr. 
Based on the CO(1-0) line, Libert et al. (2010a) derive a 
stellar systemic velocity V*.lsr = — 33.75±0.25 km s _1 , 
outflow velocities for the broad (narrow) line compo- 
nents of Vout = 17.5±0.5 km s" 1 (5.0±0.5 km s" 1 ), and 
mass loss rates for the broad (narrow) components of 
M = 6.1 x 10~ 7 M yr" 1 (1.0 x 10~ 7 M G yr" 1 ). Results 
from the CO(2-l) line were comparable to within uncer- 
tainties. We note that the outflow velocity of V1942 Sgr 
is unusually large for an Lb variable (e.g., Olofsson et 
al. 2002), a fact that may be linked with its DY Persei 
status. 

The proper motions derived from Hipparcos measure- 
ments of V1942 Sgr are 10.98±0.59 mas yr -1 in right as- 
cension and —5.10 ± 0.43 mas yr -1 in declination. From 
these values we derive a Galactic peculiar space velocity 
Vpace=40.2 km s _1 at a position angle of 95°1 (Table 5). 

4.6.1. VLA Results for V1942 Sgr 

Hi channel maps for V1942 Sgr derived from our nat- 
urally weighted VLA data are presented in Figure [T51 
An Hi total intensity map derived from the same data 
is shown in Figure Q1|J As seen in Figure [TBI ( and in Fig- 
ure[20l discussed below), V1942 Sgr lies at a velocity that 
suffers very low contamination from H 1 emission along 
the line-of-sight (see also Libert et al. 2010a). We are 
able to trace circumstellar H 1 emission at or near the 
stellar position at > 3cr significance across five contigu- 
ous spectral channels (Vlsr = —35.2 to —30.1 km s _1 ). 
The emission distribution at each of these velocities ap- 
pears somewhat amorphous, but with a sharp drop-off 
in Hi surface density at its boundaries. Versions of our 
Hi data cube made using u-v tapering (not shown) do 
not reveal any emission extended beyond what is seen in 
Figure [EH 

To first order, the velocity-integrated Hi emission dis- 
tribution of V1942 Sgr (Figure IT0|) has a shell-like mor- 
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Fig. 18. — H I channel maps bracketing the LSR velocity of V1942 Sgr, derived from naturally weighted VLA data. The spatial resolution 
is ~ 89" X 49". A star symbol indicates the stellar position. Contour levels are (— 4.2[abscnt],— 3,3,4.2,6,8.4)x2.7 mjy beam -1 . The lowest 
contour is ~ 3o\ The ficld-of-view shown is comparable to the VLA primary beam. 
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Fig. 19. — Hi total intensity map of V1942 Sgr derived from naturally weighted VLA data. The spatial resolution is ~ 89" X 49". A star 
symbol indicates the stellar position and the red line indicates the direction of space motion. The map was constructed from emission span- 
ning LSR velocities —35.2 to —30.1 km s — 1 . Intensity levels are to 120 Jy beam" 1 m s" 1 and contour levels are (1.4,2,2.8,4,5.6,8) X 14.1 Jy 
beam -1 m s — 1 . To minimize noise in the map, data at a given point were blanked if they did not exceed a 2<r threshold after smoothing 
by a factor of three spatially and spectrally. 
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Fig. 20. — Spatially integrated Hi spectrum of V1942 Sgr. The 
spectrum was derived from naturally weighted VLA data (after 
correction for primary beam attenuation) by integrating the emis- 
sion in each velocity channel within a 4.'7 X 5.'5 aperture centered at 
Qj2000=19 h 19 m 06.5 s , <5j2opo = -15°54' 45'.'0. The vertical line 
indicates the stellar systemic velocity derived from CO observa- 
tions. The dotted line shows a Gaussian fit. 

phology, reminiscent of Y UMa (JO|) and Y CVn (i}4~5j). 
Similar to those other two stars, the outer boundaries 
of the shell are irregular. The total extent of the Hi 
emission of V1942 Sgr is roughly A'.O (^0.6 pc) east- west 
and 4f6 (^0.7 pc) north-south, as measured through the 
position of the star. The stellar position is offset ~ 30" 
to the east of the geometric center of the Hi nebula. 
Other noteworthy features are the two enhancements in 
Hi surface brightness that are observed southeast and 
northwest of the star, respectively (see also Figure [TH]). 
The northwestern feature is the brighter of the two and 
lies atop a spatially extended ridge of enhanced surface 
brightness. Both of the H I enhancements lie along a line 
that bisects the stellar position at PA« 126°, although 
the significance of this is unclear, and the two surface 
brightness enhancements do not show any obvious rela- 
tionship to the position angle of the peculiar space mo- 
tion of the star (PA« 95°). 

A Gaussian fit to the integrated Hi spectrum in Fig- 
ure [20] yields a peak flux density of 0.240±0.009 Jy, a 
FWHM line width of 3.23±0.13 km s"\ and a line cen- 
troid of V LS r = -33.04 ± 0.06 km s" 1 . From these 
fit parameters, we derive an integrated Hi flux density 
/ Shi dv = 0.82±0.04 Jy km s _1 , translating to a total 
Hi mass Mhi ~ 0.055 Mq. Our measurements are in 
good agreement with the H I line parameters reported by 
Libert et al. (2010a). 

4.6.2. Discussion of V1942 Sgr Results 

The maximum angular extent of the H i-emitting shell 
that we observe surrounding V1942 Sgr (~4f6 across) 
is markedly smaller than the FIR diameter derived by 
Young et al. (1993a) from IRAS data (~6f4). Further- 
more, the good agreement between the Hi size and Hi 
line parameters that we derive in the present study com- 
pared with those derived from the NRT by Libert et al. 
(2010a) also suggests that V1942 Sgr has little or no cir- 
cumstellar H I emission extended beyond what is detected 
in our VLA images. Thus in contrast to the carbon star 
Y CVn (§ 14. 5p , our findings suggest that the Hi shell 



of V1942 Sgr may be significantly more compact than 
the FIR-emitting shell. Simulations by van Marie et al. 
(2011) have shown that this type of segregation can occur 
preferentially for large dust grains (radius <0.105/im). 
Higher resolution observations would be of considerable 
interest for confirming the FIR extent of V1942 Sgr and 
for examining in detail the relationship between the FIR 
and H I emission. 

Libert et al. (2010a) found evidence in their integrated 
Hi line profile of V1942 Sgr for a weak (~ 6±2 mJy) 
pedestal of emission extending from —39 to —27 km s — 1 . 
Based on the results of numerical modeling, they at- 
tribute this weaker, broader line component to a counter- 
part to the freely expanding wind at small radii (as traced 
by the narrower of the two CO line components detected 
by these authors). They interpret the dominant narrow 
H I line component as arising from a detached shell that 
has been decelerated as a result of its interaction with 
the ISM. Thus the interpretation of the line profile is 
analogous to the cases of RX Lep, Y UMa, and Y CVn 
(see Appendix A). 

We do not find any evidence for a counterpart to the 
broad H I emission component in our integrated VLA line 
profile (Figure [20|) or in our H I channel maps fFigure[T8|). 
In the channel maps, we would expect the broad compo- 
nent to give rise to emission at or near the stellar posi- 
tion at velocities of roughly T4 ± 5 km s _1 . However, the 
absence of this component in our data is not surprising. 
Since its expected peak line strength of only ~ 6±2 mJy, 
it is near the limits of detectability in our present data 
(where RMS noise per channel is ~2.7 mJy beam -1 ; Ta- 
ble 4). 

We have also searched our Hi data cube for a coun- 
terpart to the broad CO line component detected by 
Libert et al. (2010a). Consistent with Libert et al., 
we find no evidence of H I emission at velocities near 
Kt±17.5 km s _1 . Libert et al. suggested that this high- 
velocity outflow is likely to result from very recent mass- 
loss ( <^10 4 years) and that it may be a signpost of 
the star being close to the end of its life on the AGB. 
This possibility is interesting in light of the newly estab- 
lished DY Persei classification of V1942 Sgr, since some 
workers have suggested that these may be highly evolved 
stars nearing the planetary nebula stage (see Miller et al. 
2012 for discussion). Hi and CO observations of addi- 
tional DY Persei stars would be valuable for establishing 
whether the CSE properties of V1942 Sgr are typical of 
its class. 

We have estimated an age for the Hi shell surround- 
ing V1942 Sgr using a simple numerical model (see Ap- 
pendix A and Table Al for details). The resulting fit to 
the global H I spectrum is presented in Figure l2"gq . We 
estimate that a duration of ~ 7.2 x 10 5 yr was required to 
form the observed shell. This is ^20% higher than the 
previous estimate from Libert et al. (2010a) based on 
constraints from single dish data and more than 6 times 
larger than the CSE age derived by Young et al. (1993b) 
based on FIR data. 

4.7. AFGL 3099 

AFGL 3099 (=CRL 3099=IZ Peg) is a carbon star and 
Mira-type variable with a period of 484 days (Le Bertre 
1992). This star has a very low effective temperature 
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(1800-2000 K; Gerard & Le Bertre 2006) and is the 
coolest AGB star observed in H I with the VLA to date. 
We adopt a distance to the star of 1500 pc based on the 
period-luminosity relation of Groenewegen & Whitelock 
(1996). 

AFGL 3099 is thought to be in a highly evolved evo- 
lutionary state (Gehrz et al. 1978). It is similar to the 
prototypical carbon star IRC+10216 in that one of the 
few known examples of Galactic carbon stars to show 
regular, large-amplitude variability and exhibit a pulsa- 
tion period longer than normal (optically selected) car- 
bon Miras (Feast et al. 1985). 

Knapp & Morris (1985) detected CO(1-0) emission 
from AFGL 3099 from which they derived a rather high 
mass-loss rate of 1.3xl0 -5 Mq yr _1 (after scaling to our 
adopted distance) . They also found a stellar systemic ve- 
locity of V*,l,sr = 46-6 ± 0.4 km s _1 and a wind outflow 
speed of V out = 10.1±0.5 km s _1 . The extremely red 
color of AFGL 3099 implies it is surrounded by a dense 
circumstellar dust shell (Gehrz et al. 1978). AFGL 3099 
was not found to be an extended FIR source based on 
IRAS observations (Young et al. 1993a), although this 
may have been a consequence of the relatively large stel- 
lar distance. 

4.7.1. VLA Results for AFGL 3099 

H I channel maps of the AFGL 3099 field are presented 
in Figure [21] and a global H I spectrum coincident with 
the stellar position is presented in Figure [5H We find 
no evidence for Hi emission associated with the CSE of 
AFGL 3099 despite the presence of negligible interstellar 
contamination. No emission is seen coincident with the 
stellar position at > 2.5er in any spectral channel within 
the band. Assuming a Gaussian line profile with FWHM 
equal to twice the outflow velocity measured in CO, we 
derive a 3cr upper limit for the Hi mass of AFGL 3099 
within a single synthesized beam centered on the star of 
Mhi < 0.058M Q . This is consistent with the la limit 
previously reported by Gerard & Le Bertre (2006). 

4.7.2. Discussion of AFGL 3099 Results 

The absence of detectable H I emission directly toward 
the position of AFGL 3099 despite its high mass loss 
rate is not entirely surprising. AFGL 3099 is signifi- 
cantly more distant than other stars in our sample (Ta- 
ble 1), and the extremely low effective temperature of 
the star suggests that the bulk of its wind may be molec- 
ular (Glassgold & Huggins 1983) or simply too cold to 
emit efficiently in the 21-cm line. However, another sim- 
ilarly cool and distant carbon Mira, AFGL 3068, was 
detected in Hi by Gerard & Le Bertre (2006). Based on 
the data of Knapp & Morris (1985), AFGL 3099 has a 
~3 times weaker integrated CO line intensity compared 
with AFGL 3068, but its peak Hi intensity is more than 
8 times weaker. 

In the case of the cool and highly evolved carbon Mira 
IRC+10216 (T cff ~2200 K), Matthews & Reid (2007) 
found no H i emission coincident with the position of the 
star (see also Le Bertre & Gerard 2001), but detected arcs 
of H I emission displaced from the central position. These 
arcs were later found to be associated with the interface 
between the CSE and the ISM and to lie along an as- 
trosheath and vortical tail visible in the FUV (Matthews 



et al. 2011a). We do not find any sign of analogous fea- 
tures in the AFGL 3099 field, although the large distance 
of this star from the Galactic plane suggests that the lo- 
cal ISM density is likely to be significantly lower than in 
the case of IRC+10216. 

Another possible explanation for the absence of ex- 
tended H I emission around AFGL 3099 is that the star 
has only recently began undergoing mass loss at a high 
rate (cf. Libert et al. 2010a). We further note that the 
outflow velocity for AFGL 3099 is surprisingly low for 
a star this red (e.g., Loup et al. 1993, their Figure 8), 
suggesting that its mass-loss properties may be atypical. 

4.8. TX Piscium (TX Psc) 

The N-type carbon star TX Piscium (TX Psc) is fre- 
quently classified as an irregular (Lb) variable, although 
evidence for a period of ^220 days has been reported 
by Wasatonic (1995, 1997). TX Psc was among the first 
evolved stars in which technetium was detected (Merrill 
1956) and it is thought to be in a post-Mira evolutionary 
phase (Judge & Stencel 1991). TX Psc is the warmest 
carbon star that has been imaged in H I to date, with 
T off = 3115+130 K (Bergeat et al. 2001). Based on the 
Hipparcos parallax of 3.63+0.39 mas (van Leeuwen 2007) 
we adopt a distance to TX Psc of 275 pc. A recent anal- 
ysis by Klotz et al. (2013) suggests that the initial mass 
of TX Psc was in the range 1-3 Mq. 

Using the CO(1-0) data of Eriksson et al. (1986), Loup 
et al. (1993) derived an outflow velocity for TX Psc's 
wind of V ou t = 10.7 ± 1.9 km s" 1 and a mass loss rate 
M = (3.1 ± 1.6) x 10~ 7 Mq yr^ 1 (scaled to our adopted 
distance). The stellar systemic velocity as derived from 
CO observations varies among different authors depend- 
ing on the telescope and the transition used; here we 
adopt the mean of the values quoted by Olofsson et al. 
(1993): V* jL sr = 12.2+1.1 km s -1 , where the quoted un- 
certainty is based on the dispersion in the measurements. 
The Hipparcos proper motions of —33.68 ± 0.40 mas 
yr _1 in right ascension and —24.49 ± 0.32 mas yr _1 in 
declination imply a peculiar space velocity of V^pacc = 
66.4 km s _1 along a position angle of 248° (Table 5). 

4.8.1. VLA Results for TX Psc 

We present H I channel maps for TX Psc in Figure [23J 
In the spectral channel centered at Vlsr =12.4 km s _1 
(which is the nearest to the systemic velocity of the star), 
H I emission is detected at the position of TX Psc and is 
seen extending several arcminutes to the northeast. Spa- 
tially resolved emission is also clearly detected toward 
the stellar position at Vlsr = 13.6 km s _1 . The spectral 
channels with Vlsr > 12-4 km s _1 appear to suffer only 
minimally from line-of-sight confusion, implying that the 
detected emission almost certainly arises from the CSE 
of TX Psc. The distribution of this material is consis- 
tent with a gaseous wake trailing the star. The H I emis- 
sion structure in these velocity channels also resembles 
the morphology and orientation of the faint FIR-emitting 
wake seen in the Herschel images of Jorissen et al. (2011) 
and the Spitzer 70/im image of Geise (2011). 

Looking toward more blueshifted velocities (Vlsr =7.2 
to 11.1 km s _1 ), we see evidence of an increasingly ex- 
tended trail of emission stretching behind TX Psc. This 
emission appears well-collimated out to roughly 8' from 
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Fig. 21. — Hi channel maps bracketing the systemic velocity of AFGL 3099 derived from naturally weighted VLA data. The spatial 
resolution is ~ 64" X 52". A star symbol indicates the stellar position. Contour levels are (— 3,3,4. 2[absent]) X 1.7 mjy beam -1 . The lowest 
contour is ~ 3o\ The field-of-view shown is comparable to the VLA primary beam. 



the star. However, within this velocity range there are 
also clear signatures of large-scale emission that is poorly 
sampled by the VLA. The incomplete spatial sampling of 
this large-scale emission results in the occurrence of posi- 
tive and negative bands across the field (see also §4.8.2p . 
These artifacts make it difficult to accurately and un- 
ambiguously determine the full extent of the H I wake 
associated with TX Psc. However, a careful examination 
of the data does provide some important clues. 

In the spectral channel centered at Vlsr = 
11.1 km s^ 1 , a "kink" is visible in the emission trailing 
to the northeast of TX Psc, near a J20 oo = 23 h 46 m 48 s , 
<5j20oo = 03°33'02". The presence of this discontinuity 
raises some doubt that the material east of this position 



is related to the circumstellar wake. Furthermore, in the 
Vlsr =8.5 and 9.8 km s _1 velocity channels, an extended 
H I "cloud" is observed to the northeast of the star whose 
characteristics suggest that it is not associated with the 
CSE of TX Psc. This becomes further evident in the H I 
total intensity map in Figure l24l 

Figure |2"41 highlights a tapering of the emission trailing 
TX Psc, occurring roughly 8' behind the star. This is 
analogous to what has been observed in the H I wakes of 
other AGB stars (Matthews et al. 2008, 2011b). How- 
ever, beyond this, there is a dip in the Hi column density 
followed by an abrupt change in the opening angle of the 
trailing emission. The marked discontinuity in column 
density is evident even in the u-v tapered version of our 
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Fig. 22.— VLA Hi spectrum of AFGL 3099 obtained by aver- 
aging over a single synthesized beam centered on the star. The 
vertical line indicates the systemic velocity of the star as deter- 
mined from CO observations. 

data where we expect improved sensitivity to extended, 
low surface brightness emission. Furthermore, the over- 
all morphology of this emission beyond this point (i.e., at 
r >8' from the star) appears to be generally inconsistent 
with a circumstellar wake (cf. Villaver et al. 2002, 2012; 
Wareingetal. 2007a, b, c). While wide-angle circumstel- 
lar wakes may indeed arise from mass-losing stars under 
certain conditions (e.g., Comeron & Kaper 1998; Ware- 
ing et al. 2007a; Sahai & Chronopoulos 2010; Matthews 
et al. 2012), a change from narrow collimation near the 
star to a broad, wide-angle wake at larger projected dis- 
tances would be difficult to explain. For these reasons, 
we conclude that only the gas within <8' from TX Psc 
can be convincingly associated with a trailing circum- 
stellar wake. While we cannot rule out the presence of 
some additional circumstellar debris beyond this radius 
(see below), the dominant emission in our VLA maps at 
these projected distances from the star appears to be in- 
terstellar in origin. Under these assumptions, we have 
derived the spatially integrated H I spectrum of TX Psc 
shown in (Figure l25j). Based on a Gaussian fit to this 
spectrum (including only data with Vlsr > 8.5 km s _1 ), 
we derive a peak Hi flux density of 0.35±0.02 Jy, a line 
centroid of Vlsr = H.4±0.1 km s -1 , and a FWHM line 
width of 4.0 ±0.2 km s _1 . From this we measure an in- 
tegrated Hi flux density of J Shi dv = 1.48 ± 0.11 Jy 
km s" 1 , translating to an Hi mass A/hi ~ 0.026 Mq. 

The wake of emission trailing behind TX Psc lies along 
a position angle of ~ 45°, and is thus displaced by —23° 
compared with a bisector defined by the space motion 
vector of the star (see Table 5). This suggests the pres- 
ence of a transverse flow in the local ISM (see below). 

Perhaps the most striking feature of the H I emission 
toward TX Psc is its velocity field (Figure [26l) . While 
the velocity field map contains a great deal of spurious 
structure resulting from undersampled large-scale emis- 
sion, some noteworthy features are nonetheless apparent, 
including velocity gradients both parallel and perpendic- 
ular to the space motion of the star. These are discussed 
in detail below. 

4.8.2. Discussion of TX Psc Results 



TX Psc was previously observed in Hi by Gerard & 
Le Bertre (2006). Its CSE was clearly detected, although 
the blue edge of the line profile suffered from confusion. 
Gerard & Le Bertre found evidence that the H I emission 
associated with TX Psc is quite extended, with a radius 
of ~ 12'. This is significantly larger than the 60/mi ra- 
dius of 3!l reported by Young et al. (1993a) based on 
IRAS data. Similarly, based on higher resolution Her- 
schel 70/Ltm data, Jorissen et al. (2011) found a source 
with a total extent of roughly 3' (as estimated from their 
Figure 6), which includes a broad, low surface brightness 
trail of emission extending toward the northeast. Geise 
(2011) measured a slightly larger angular extent (~4!6) 
from a Spitzer 70^m image. 

Examination of H I survey spectra from Kalberla et al. 
(2003) toward the direction of TX Psc reveals that there 
is a line component with peak brightness temperature 
~2 K centered near Vlsr ~10 km s _1 , very near the 
systemic velocity of TX Psc. This emission is too strong 
to be plausibly associated with the CSE of TX Psc, and 
instead is likely to be interstellar in origin. This emis- 
sion reveals its presence in our VLA maps as the pos- 
itive and negative banding seen in our channel images 
from Vlsr = 8.5 to 11.1 km s _1 (Figure [23]) implying 
that it is extended over scales of >15' (the maximum 
angular scales that can be sampled with the VLA in D 
configuration). The presence of this additional emission 
component so close in projected position and velocity to 
TX Psc suggests that the amount of extended emission in 
the CSE may have been overestimated from the NRT ob- 
servations. Indeed, the Hi mass for the CSE of TX Psc 
derived from the NRT data is nearly five times larger 
than we derive from the VLA. While we cannot rule out 
that the H I wake of TX Psc may include a low surface 
brightness extension that cannot be readily detected with 
the VLA (analogous to the one confirmed by the NRT 
toward Mira; see Matthews et al. 2008), our VLA maps 
suggest that confusion may preclude a definitive test of 
this possibility. 

A systemic decrease in the radial velocity of gas in ex- 
tended circumstellar wakes parallel to the direction of 
motion is expected as a consequence of the gradual de- 
celeration of the wake material owing to its interaction 
with the ISM (e.g., Wareing et al. 2007c; Raga & Canto 
2008). Indeed, such gradients have been detected in the 
extended wakes of Mira (Matthews et al. 2008) and 
X Her (Matthews et al. 2011b). In the case of TX Psc, 
we likewise see evidence for a systematic blueshifting of 
the "tail" gas with increasing distance from the star in 
the H I channel maps (Figure |2"3")) and H I velocity field 
(Figure [26]). The gradient is maximized along a posi- 
tion angle of ~ 58° and has a magnitude of ~l-2 km s _1 
projected in the plane of the sky. 

While in principle, quantifying the velocity gradient 
along a gaseous wake provides a means of age-dating the 
stellar mass loss history (Matthews et al. 2008, 2011b; 
Raga & Canto 2008), our ability to do this for TX Psc 
is compromised by the presence of a second prominent 
velocity gradient perpendicular to the direction of space 
motion of the star. The magnitude of this gradient (pro- 
jected in the plane of the sky) is ^4.5 km s _1 across the 
width of the Hi envelope (~ 6' or 0.48 pc). This is fur- 
ther illustrated in Figure 1271 where we show a position- 
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Fig. 23. — Hi channel maps bracketing the systemic velocity of TX Psc derived from naturally weighted VLA data. The spatial resolution 
is ~ 72" X 49". A star symbol indicates the stellar position. Contour levels are (-6, -4.2, -3, 3, 4.2, ...12) x2 mjy beam" 1 . The lowest 
contour is ~ 3cr. The ficld-of-view shown is comparable to the VLA primary beam. 



velocity cut across the emission. This perpendicular ve- 
locity gradient appears to be real and not an artifact 
caused by unresolved large-scale structure or its sidelobes 
in our maps. The most compelling evidence of this is that 
the velocity gradient persists even at velocities that are 
uncontaminated by confusion (Vlsr > 12 km s _1 ). It 
also persists if we impose u-v restrictions on the data 
during imaging in order to suppress the short spacings 
(u, v < 350A) that give rise to the spurious large-scale 
structure in the images. Unfortunately, such u-v restric- 
tions also suppress much of the large-scale emission that 
appears to be genuinely associated with TX Psc. 

For the semi-regular variable X Her, Matthews et al. 
(2011b) also reported evidence of a velocity gradient or- 



thogonal to the direction of space motion in their Hi 
data. In the X Her case, this orthogonal axis corre- 
sponds to the symmetry axis of a bipolar outflow pre- 
viously detected at smaller radii through molecular line 
observations. Matthews et al. suggested that the H 1 may 
therefore be tracing an extension of this outflow out be- 
yond the molecular dissociation radius. Might a similar 
situation apply to the case of TX Psc? 

TX Psc has no published molecular line imaging obser- 
vations, so it is unknown whether the star has an axisym- 
metric molecular outflow. Single-dish mapping by Hcskc 
et al. (1989) revealed hints of deviations from spheri- 
cal symmetry, and the CO profile published by Hcskc 
(1990) shows a two component line structure, which is 



32 



Matthews et al. 



100 



03 45 



35 



O 

H 

< 



LU 
Q 




15 
23 4715 



46 45 



30 15 00 45 45 

RIGHT ASCENSION (J2000) 



Fig. 24. — Hi total intensity map of TX Psc, derived from tapered VLA data, incorporating emission at LSR velocities from 8.5 to 
13.6 km s — 1 . The spatial resolution is ~ 104" X 79". A star indicates the stellar position and the red line indicates the direction of space 
motion. The contour levels are (-5.6, -4, -2.8, -2, -1.4, -1, 1...5.6)x28 Jy beam" 1 m s" 1 and intensity levels are 20 to 125 Jy beam 
m s — 1 . The data have been corrected for primary beam attenuation. The emission associated with the CSE of TX Psc exhibits a head-tail 
morphology, although the bulk of the emission east of aj2000 = 23 h 46 m 48 s and north of 5j2000 = 03°35'35.8" is probably unrelated to the 
star (see text). To minimize noise in the map, data at a given point were blanked if their absolute value did not exceed a 3<r threshold 
after spatially smoothing the data by a factor of three in velocity and factor of 5 spatially. The positive and negative bands running across 
the image arise from large-scale interstellar emission that is poorly spatially sampled by the VLA. 



frequently a hallmark of stars with bipolar outflows (e.g., 
Kahane & Jura f996; Bergman et al. 2000; Josselin et 
al. 2000; Libert et al. 2010b). The CO mapping data 
of Heske et al. also show evidence for a spatial offset 
between the blue and red portions of the line emission 
along a position angle roughly parallel to the direction 
of space motion of the star. However, this velocity shift 
has the opposite sign to the gradient seen in the H I along 
this direction (i.e., with blueshifted CO emission peak- 
ing to the southwest and redshifted CO emission to the 
northeast) Q Another significant problem with this expla- 
nation for the velocity gradient transverse to the TX Psc 
wake is that the velocity gradient appears to persist well 
beyond the vicinity of the star, where the gas should be 

9 We assume that right ascension increases to the left in Figure 3 
of Heske et al. 1989. However, the signs of the offsets in the x-axis 
labeling create an ambiguity. 



dominated by a turbulent flow of ram pressure stripped 
material, not the stellar outflow. 

An alternative explanation for the velocity gradient or- 
thogonal to its H I wake is the influence of a transverse 
flow in the local ISM. Evidence for local transverse flows 
has now been seen in the neighborhood of several other 
evolved stars (e.g., Ueta et al. 2008; Ferguson & Ueta 
2010; Ueta et al. 2010; see also = j4.2.2|) . although the ori- 
gin of these flows remains unclear. In the case of TX Psc, 
such a transverse flow could also offer a natural explana- 
tion for the offset in position angle between the observed 
H I wake and the derived space motion vector in Fig- 
ure El 

We can attempt to quantify the magnitude and direc- 
tion of a local ISM flow following the procedure outlined 
in Ueta et al. (2008). Wilkin (1996) derived an analytic 
expression that relates the shape of a bow shock surface 
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Fig. 25. — Spatially integrated Hi spectrum of TX Psc derived 
from naturally weighted VLA data (after correction for primary 
beam attenuation). Emission integrated over two adjacent rectan- 
gular apertures was summed to produce the final spectrum; the first 
aperture was 6.'2x6.'8, centered at o J2 ooo = 23 h 46 m 29.8 s , <5 J2 ooo = 
03°30'27.5"; the second was 2.'8x4.'0, centered at a J20 oo=23 h 46 m 
47.9 s , <5j2000 = 03°32'42.5". The portion of the spectrum blue- 
ward of Vlsr ~& km s _1 is affected by sidelobes from large-scale 
Galactic emission. The vertical bar indicates the systemic velocity 
of the star as determined from CO observations. 



to its standoff distance, Rq: 

R(6) = i? csc 9y/3(T 



cot6>) 



(1) 



where 9 is the angular displacement from the bow shock 
apex as measured from the position of the star, and: 



MVo, 
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space 



Here mn is the mass of the hydrogen atom, njj is the 
H I number density of the local ISM, and [i is the mean 
molecular weight of the gas (assumed to be 1.3). Rq effec- 
tively defines the locus of ram pressure balance between 
the internal pressure of the CSE surrounding a mass- 
losing star and the ram pressure exerted by the ISM due 
to the star's supersonic space motion. 

Based on a three-dimensional fit of Equation 1 to the 
Herschel FIR data of TX Psc, Cox et al. (2012a) de- 
rived a deprojected standoff distance, Ro — 37" (1.52 x 
10 17 cm), with a position angle and inclination for the 
bow shock of 238° and 11°, respectively. Inverting Equa- 
tion 2 and using the fitted Rq value together with M 
and V au t from Table 1, one may then derive the peculiar 
space velocity of the star with respect to its local ISM: 
K'pacc^SS.Gn^kms- 1 . 

Adopting the position angle and inclination of the bow 
shock from Cox et al., V.' pace can now be decomposed into 
its equatorial space- velocity components: 
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From the measured proper motion of TX Psc (see above) , 
the heliocentric Galactic space velocity of the star in 
equatorial coordinates (see Table 5) can be expressed as 
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The local ISM flow in the vicinity of TX Psc may then 
be obtained as the difference between V£ pace and Vpace : 
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Based on the distance of TX Psc from the Galactic 
plane (—212 pc; Table 5), the expected ISM density in 
the vicinity of TX Psc is n H ~0.2 cm~ 3 (Table 5). Adopt- 
ing this value yields Vnow ~ 60 km s _1 along a position 
angle PA=338°. While the true value of is uncertain, 
the nominal value predicts a position angle for the flow 
that corresponds precisely to the direction of the vertical 
velocity gradient across the TX Psc envelope ( Figure [27l). 
It thus seems quite plausible that the velocity gradient 
across the TX Psc shell results from a transverse ISM 
flow rather than outflow from the star itself. The direc- 
tion of the observed velocity gradient is also consistent 
with this interpretation. 



5. GENERAL DISCUSSION 



5.1. 



A Synopsis of Hi Imaging Results to Date 

The present study now brings the total number of red 
giants successfully imaged in the Hi 21-cm line to 12. 
These comprise 11 AGB stars and the red supergiant 
a Ori (see also Bowers & Knapp 1987, 1988; Matthews 
& Reid 2007; Matthews et al. 2008, 2011a, b; Le Bertre 
et al. 2012). Including the current study, five additional 
stars have also been observed with the VLA but were not 
detected (see also Hawkins & Proctor 1993; Matthews & 
Reid 2007) E3 This is only a small fraction of the total 
sample of evolved stars now surveyed in the H I line from 
single-dish surveys. For a more extensive discussion of 
Hi detection statistics as a function of stellar properties, 
we therefore direct the reader to the single-dish surveys 
of Gerard & Le Bertre (2006) and Gerard et al. (2011b). 
Nonetheless, the imaged sample is now large enough to 
enable investigation of some emerging trends and to be- 
gin to explore how H I properties correlate with other 
properties of the stars. 

We emphasize that the current sample is a biased one 
in that imaging studies to date have largely (although 
not exclusively) targeted stars that were already securely 
detected in Hi, based on single-dish surveys (see §2; 
Matthews & Reid 2007). Additionally, since a portion 
of the Hi imaging sample was chosen to be displaced 
from the predominant Galactic H I emission along the 
line-of-sight, the sample may contain a disproportionate 
number of stars with moderately high space velocities, 
whose ages and metallicities may not be representative 
of mean Galactic values. 

10 Knapp & Bowers 1983 presented upper limits on the circum- 
stellar Hi content of 16 evolved stars based on VLA snapshot ob- 
servations; however, several of these stars have been subsequently 
detected, and the upper limits for the remainder cannot rule out 
the presence of circumstellar Hi at comparable levels. 
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Fig. 26. — Intensity-weighted Hi velocity field of TX Psc, derived from tapered VLA data. The spatial resolution is ~ 104" X 79". The 
arrow indicates the direction of space motion of the star. Note th e pr esence of velocity gradients both parallel and perpendicular to the 
direction of motion. The H I total intensity contours from Figure 1241 (positive values only) are overplotted for reference. The color bar 
shows LSR radial velocity in units of km s . 



One unifying trend among the stars successfully im- 
aged in Hi so far is that in all cases, the Hi images and 
spatially resolved kinematics display evidence for inter- 
action between the CSE and surrounding medium. This 
local medium includes the ISM, but may also include 
material shed during earlier mass-loss episodes. Mani- 
festations of this interaction come in several forms, in- 
cluding: (1) extended wakes, trailing the motions of the 
stars through space (hereafter "Category 1"); (2) quasi- 
stationary detached shells, whose properties can be mod- 
eled as the slowing of ejected wind material through its 
interaction with the ISM (hereafter "Category 2"); (3) 



arcs or clumps of gas displaced from the stellar position 
that correlate with bow shocks or other signatures of in- 
teraction with the ISM seen at other wavelengths (here- 
after "Category 3"). A few stars display features of more 
than one of these classifications. 

To date, five stars have been imaged that display 
the trailing Hi wakes characteristic of Category 1: 
RX Lep, TX Psc, X Her, Mira, and RS Cnc "(see 
gj2 03] Matthews & Reid 2007; Matthews et al. 
2008, 2011b; Libert et al. 2010b). Among these stars, 
four have in common relatively high space velocities 
(^space >57 km s _1 ). This is consistent with the results 
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TABLE 5 
Derived Properties of the Sample 



Name 


^pcak 


JSdv 


Vlsr,hi AVhi 


M m 


8h 


Vspacc 


PA 




"HI 




(Jy) 


(Jy km s ) 


(km s _1 ) (km s — *) 


(Mq) 


(') 


(km s — 1 ) 


(°) 


(PC) 


(cm -3 ) 


(1) 


(2) 


(3) 


(4) (5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 




Oxygen-rich stars 


IK Tau* 


<0.006 


<0.10 




< 0.001 




>34.6 




-113 


0.65 


RX Lep 


0.58±0.02 


2.93±0.22 


28.4±0.1 5.0±0.3 


0.015 


17.'5 


56.6 


24.6 


54 


1.16 


Y UMa 


0.10±0.01,0.03±0.01 


0.63±0.13 


16.4±0.2,17.2±0.8 3.2±0.4,9.2±2.3 


0.022 


6.'8 


19.2 


57.1 


352 


0.06 


R Peg 


0.16±0.01 


0.76±0.09 


22.4±0.2 4.4±0.4 


0.029 


£15' 


26.3 


147.6 


-266 


0.14 


Carbon-rich stars 


Y CVn 


0.79±0.02 


2.68±0.11 


20.65±0.05 3.2±0.1 


0.047 


9' 


31.0 


36.8 


273 


0.13 


V1942 Sgr 


0.240±0.009 


0.82±0.04 


-33.04 ±0.06 3.2±0.1 


0.055 


4.'6 


40.2 


95.1 


-107 


0.69 


AFGL3099* 


<0.0024 


<0.11 




< 0.058 




>46.6 




-1081 


0.00 


TX Psc 


0.35±0.01 


1.48±0.11 


11.4±0.1 4.0±0.2 


0.026 


~8' 


66.4 


248.0 


-212 


0.24 


Note. — Explar 


ation of columns: (1) target n£ 


me; (2) peak H I 


flux density in the spatially integrated line profi 


c; upper limits 


ire 3r/; 


(3) integrated 


flux density 


of the Hi 1 


nc profile, 


based on a Gauss 


an fit; (4) LSR velocity of the 


Hi line center be 


sed on a Gaussian fit; two values are quoted for 


multi-componcr 


t lines 


(5) FWHM lin 


ewidth of th 


e H I profil 


e based on 


a Gaussian fit; two values arc quoted for multi-t 


omponent lines; 


(6) Hi mass of the CSE derived from the currcr 


t VLA observat 


ons; u 


3pcr limits arc 


3<t; (7) «» 


imum angular extent 


of the observed ci 


rcumstollar Hi emission; (8) p 


=culiar spaec vel, 


city (in LSR frame); (9) position angle of space 


motion project 


d on the plane of the 


sky (measu 


red east fr 


Dm north), 


derived using the 


prescription of Johnson & Sod 


erblom 1987 and 


the solar constants from Schonrich et al. 2010; 


(10) distance from the 


Galactic plane 


, estimated 


as z dsir 


b + 15 pc, 


where values of d 


stance d and Galactic latitud* 


6 are taken fro 


■n Table 1; (11) local ISM number density, csim 


ated from n-^(z 


) = 2.0e lODpc ( Lo 


up et al. 1993), where 




from column 10. 



















Proper motion measurements arc unavailable; no space motion vector was derived. 



five stars have mass loss rates within a narrow range of 
values M w (1 - 3) x 10~ 7 A/ Q yr _1 , despite a range of 
outflow velocities and variability classes. 

The current study has identified four stars belonging 
to Category 2 (i.e., that appear to be surrounded by 
quasi-stationary detached shells): RX Lep ( fc|4.2[) . Y UMa 
(MM, Y CVn (S3]), and V1942 Sgr ( U^ . The red su- 
pergiant a Ori (Le Bertre et al. 2012) also falls into 
this category. All four of the AGB stars in this cate- 
gory have mass-loss rates in the same narrow range as 
the stars with Hi tails [M « (1 - 3) x 10~ 7 M© yr -1 ]. 
However, the space velocities of these stars (Table 5) are 
on average lower, ranging from 19.2 km s^ 1 for Y UMa 
to 56.6 km s _1 for RX Lep (which also exhibits an ex- 
tended Hi tail). Similarly, a Ori has T4 pacc ~31 km s _1 
(Le Bertre et al. 2012). Our imaging confirms that 
detached Hi shells with similar properties can occur in 
both carbon-rich and oxygen-rich stars. As already sug- 
gested by Libert et al. (2007), the formation mechanism 
for these H I shells is almost certainly different from the 
thin molecular shells seen around some carbon stars (cf. 
Schoier et al. 2005; Olofsson et al. 2010). Some key dif- 
ferences are that the molecular shells are extremely thin 
and nearly perfectly circular, whereas the Hi shells are 
geometrically thick and have shapes ranging from irreg- 
ular to boxy to egg-shaped. 

We have presented the results of simple numerical mod- 
eling of the four AGB H I shells from which we estimate 
ages in the range 9.0 x 10 4 to 7.2 x 10 5 yr (Appendix A; 
see also Le Bertre et al. 2012 for the case of a Ori). 
For Y UMa, the age that we derive for the CSE is com- 
parable to that previously derived from FIR observa- 
tions by Young et al. (1993b). However, for Y CVn 
and V1942 Sgr (both carbon stars) the duration of the 
mass loss history inferred from our new H I observations 
is nearly 7 times larger than derived previously from FIR 
observations. Indeed, as more H I measurements of the 
extended CSEs become available, they continue to point 
toward a past systematic underestimation of the mass 
loss timescales for many red giants (see also Matthews et 
al. 2008, 2011b; Le Bertre et al. 2012). 
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Fig. 27. — Position-velocity cut perpendicular to the direction 
of spa ce m otion of TX Psc, derived from velocity field shown in 
Figure 1261 The cut was extracted along a position angle of 158° 
and is centered northeast of the stellar position, at aj2000=23 h 
46 m 26.8 s , <5j2000 = 03°31'17.5". Negative displacements are to 
the southeast and positive displacements to the northwest. 



of hydrodynamic simulations (e.g., Wareing et al. 2007a; 
Villaver et al. 2012) that show that high space velocity 
is one of the primary criteria leading to the development 
of extended gaseous tails. The one exception is RS Cnc, 
with V spacc wl2.4 km s" 1 . This star is located relatively 
close to the Galactic plane (z «97 pc), so its Hi tail 
may be enhanced through the effects of a high local ISM 
density. Simulations also show that cometary tails may 
form very early in the AGB stage even for stars with low 
space velocities (Villaver et al. 2012). Furthermore, since 
the proper motion values for this star have large uncer- 
tainties (van Leeuwen 2007), the space velocity may be 
underestimated. 

One feature in common to all of the stars so far found 
to have extended gaseous wakes is that their mean effec- 
tive temperatures are high enough (T c g >2500 K) such 
that their winds are predicted to be dominated by atomic 
hydrogen (Glassgold & Huggins 1983). In addition, all 
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One of the stars recently identified as belonging to 
Category 3 (i.e., stars surrounded by arcs or clumps 
of gas displaced from the stellar position) is a Ori 
(Le Bertre et al. 2012), while the other two stars we 
place in this class, IRC+10216 and R Cas (Matthews 
& Reid 2007; Matthews et al. 2011a), are Mira vari- 
ables. All three of these stars have high mass-loss rates 
(M > 1 x 10~ 6 M Q yr" 1 ), although a Ori has a much 
higher effective temperature (T c g = 3641 ± 53 K; Perrin 
et al. 2004) compared with the two Miras (T cS <2500 K). 
For both IRC+10216 and R Cas, the Hi imaged by the 
VLA comprises isolated arcs of emission displaced from 
the stellar position, although the two stars differ dramat- 
ically in the size scales of their CSEs (cf. Matthews & 
Reid 2007), and they are believed to be at very differ- 
ent evolutionary stages. In the case of IRC+10216, it 
is suspected that the observed H I emission is partly the 
result of swept-up interstellar material, although there 
is evidence that it is augmented by photodissociated H2 
(Matthews et al. 2011a). 

Only one star imaged in H 1 to date, the Mira variable 
R Peg, does not fall neatly into any of the three classifi- 
cations defined above. While its H 1 mass and integrated 
Hi profile are similar to the other detected stars, its mor- 
phology is rather peculiar ( ^4.4|) . R Peg has a higher 
mass loss rate than any of the other VLA sample stars 
where Hi has been detected directly toward the stellar 
position (M «5x 10~ 7 M q yr -1 ), suggesting that it may 
be somewhat more evolved. In this case, one possibility 
is that the ejecta that we observe in H I are left over from 
previous cycles of mass loss and have been subsequently 
distorted by transverse flows in the surrounding ISM. 

For all stars detected in Hi so far, we have esti- 
mated the number density of hydrogen atoms in the lo- 
cal ISM using the expression from Loup et al. (1993): 

n-ft(z) = 2.0e~ 10 °p c cm -3 , where z is the distance from 
the Galactic Plane in pc, estimated as z ~ cfein(6) + 15 pc. 
Here d is the distance to the star and b is the Galactic 
latitude. For the full VLA-detected sample, the values 
range from nn ~ 0.059 cm -3 for Y UMa (Table 5) to 
7iH ~ 1-7 cm -3 for a Ori and R Cas. For the Category 1 
stars, values of njj range from 0.24 cm -3 for TX Psc to 
1.2 cm -3 for RX Lep (Table 5). Thus there is no evi- 
dence that trailing H 1 wakes are preferentially found as- 
sociated with stars in dense local enviroments, although 
as noted in i )4.2.2[ the moderately high local ISM density 
surrounding RX Lep may help to account for its highly 
collimated wake. We also emphasize that the nn values 
computed using the above formula are rather uncertain 
owing to their dependence on the adopted Galactic model 
and small-scale fluctuations in the ISM density across the 
Galaxy. 

The sample of five stars so far undetected in H 1 with 
the VLA is still too small to draw any firm conclu- 
sions, but here we briefly review their properties to help 
guide future studies. First, we point out that the two 
undetected stars from the current study, IK Tau and 
AFGL 3099, have temperatures and mass loss rates com- 
parable to two previously detected Category 3 Miras: 
R Cas and IRC+10216. However, Gerard et al. (2011b) 
have found that in general, cooler stars with higher mass 
loss rates are less frequently detected in Hi. This may 



be due to a combination of a predominantly molecular 
composition in the CSEs, their geometrically large sizes, 
and/or temperatures in the CSE that are too cool to 
allow the detection of H 1 in emission (see £ 14. 1 .2[) . 

Having a low mass loss rate does not necessarily imply 
that a star should be more difficult to detect in Hi, since 
available M estimates generally represent only recent 
mass loss and are derived from tracers at much smaller 
radii than the extended circumstellar material traced by 
Hi. However, to date, the two stars observed with the 
VLA having the lowest mass loss rates both remain un- 
detected. The Mira R Aqr (M = 6 x 10~ 8 Mq yr" 1 ; 
Spergel et al. 1983) is part of a symbiotic binary and 
is surrounded by a hot ionized region. Matthews & 
Reid (2007) placed an upper limit on its Hi mass of 
Mm < 4.9 x 10" 4 Mq. The undetected semi-regular vari- 
able W Hya (Hawkins & Proctor 1993) also has a com- 
parably small mass loss rate (M ~ 7 x 10 -8 Mq yr -1 ; 
Olofsson et al. 2002). In this case, the lack of detec- 
tion might be related to the moderately cool stellar tem- 
perature (T eff «2450 K; Feast 1996) or the difficulty of 
disentangling weak and clumpy CSE emission from line- 
of-sight contamination (see Hawkins & Proctor 1993). 

The remaining star so far undetected in H 1 by the 
VLA is the semi-regular variable EP Aqr, whose mass 
loss rate is an order of magnitude higher than R Aqr or 
W Hya (M w 6 x 10~ 7 Mq yr" 1 ; Winters et al. 2007). 
It poses an interesting case, since it is not a known bi- 
nary, and its effective temperature is high enough to 
predict the presence of atomic hydrogen in its atmo- 
sphere (T cf f w 3240 K; Dumm & Schild 1998). The 
non-detection of this star with the VLA is particularly 
puzzling, as Le Bertre & Gerard (2004) reported an Hi 
detection of EP Aqr with the NRT, and the star is also 
an extended FIR source (Young et al. 1993a; Cox et 
al. 2012a). However, Matthews & Reid (2007) found 
only a few isolated clumps of Hi in the region, none 
of which could be unambiguously linked with the CSE. 
Deeper observations of this star with improved u-v cov- 
erage might be enlightening, as some of the clumps de- 
tected by Matthews & Reid might turn out to be density 
enhancements in an extended, diffuse envelope or pock- 
ets of photodissociated gas within an ancient molecular 
envelope. 

5.2. Future Prospects 

It is clear that Hi imaging studies of larger samples 
of stars are needed to enable better statistical compar- 
isons between the properties of the stars, their CSEs, and 
their interstellar neighborhoods. However, despite the 
small size of the current VLA H 1 sample, certain trends 
are already clear. The results to date underscore initial 
findings from single-dish surveys and unambiguously es- 
tablish the importance of the interstellar environment on 
the large-scale properties of CSEs. While additional ex- 
amples of stars interacting with their environments have 
also been identified in recent FIR and FUV surveys (e.g., 
Martin et al. 2007; Stencel 2009; Sahai & Chronopoulos 
2010; Ucta 2011; Cox et al. 2012a), Hi data have the 
unique advantage that they supply kinematic informa- 
tion on the extended circumstellar material. As demon- 
strated in the present study, such data in turn permit 
numerical modeling of the interaction between circum- 
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stellar shells and their surroundings and the determina- 
tions of the duration of the stellar mass loss history (see 
also Matthews et al. 2008, 2011b; Libert et al. 2007, 
2010a, b). Hi data can also provide crucial morpho- 
logical and kinematic information for testing hydrody- 
namic simulations of individual mass-losing stars moving 
through the ISM — a potential that has largely gone un- 
tapped (cf. Villaveret al. 2012; Wareing 2012; Mohamed 
et al. 2012). 

Another important lesson that has emerged from Hi 
imaging surveys to date is that no single tracer can pro- 
vide a complete picture of the complex circumstellar en- 
vironments of mass-losing evolved stars. If one hopes to 
understand the complete stellar mass loss history, it is 
critical to study the CSE on both large and small spatial 
scales. Moreover, it has become clear that even when 
different tracers span overlapping spatial scales (e.g., Hi 
and FIR-cmitting dust), each supplies unique yet comple- 
mentary information that can be used together to provide 
a more comprehensive picture of the kinematics, physical 
conditions, and chemistry across the entire CSE. 

6. SUMMARY 

We have presented H I 21-cm line imaging observations 
of a sample of eight AGB stars obtained with the VLA. 
We have unambiguously detected H I emission associated 
with the CSEs of six of the stars. All six of these stars 
were already known to have atomic hydrogen in their 
CSEs based on the results of single-dish surveys. How- 
ever, our imaging results reveal a wealth of new informa- 
tion on their spatially resolved morphologies and kine- 
matics. In all cases we identify clear evidence that the 
properties of the extended CSEs are significantly affected 
by interaction with their local interstellar environments. 
The detected stars have circumstellar Hi masses in the 
range Mjji ~ 0.015 — 0.055M Q and the size scales of their 
circumstellar Hi emission range from ^0.6 to 1.7 pc. 

Four stars in the present study (RX Lep, Y UMa, 
Y CVn, and V1942 Sgr) are found to be surrounded by 
quasi-stationary detached hydrogen shells whose origin 
can be explained as a result of the stellar outflow be- 
ing abruptly decelerated at a termination shock that oc- 
curs at the interface between the circumstellar and inter- 
stellar environments. These shells are seen around both 
oxygen-rich and carbon-rich stars, and appear to have a 
different origin from the thin detached molecular shells 
seen around some carbon stars. All four Hi shell stars 
have similar mass-loss rates [M ss (1 — 3) x 1O~ 7 M yr _1 ] 
and low-to-moderate space velocities, but otherwise show 
a range in properties, including chemistries, wind out- 
flow velocities, and variability classes. All of the shells 
have a clumpy rather than smooth structure, and none 
of the observed H I shells are spherical; they range from 
egg-shaped (RX Lep) to boxy (Y CVn) to amorphous 
(Y UMa and V1942 Sgr). In the case of Y UMa, there is 
evidence that one edge of the shell was distorted through 
its interaction with the local ISM. The shell of RX Lep 
exhibits a pronounced elongation along the direction of 
space motion of the star, as well as a "cracked egg" mor- 
phology, with a dearth of emission along its equatorial 
region. 

Based on simple numerical models we have estimated 
ages for the four shells in the range 9.0 x 10 4 to 7.2 xlO 5 yr 
(Appendix A). These results support previous sugges- 



tions that RX Lep is in a relative early evolutionary state 
compared with the other three stars. For the carbon stars 
Y CVn and V1942 Sgr the derived ages also imply that 
mass loss time scales previously derived from FIR mea- 
surements have been underestimated by a factor of ~7. 

For both Y UMa and Y CVn, high-resolution FIR im- 
ages are available that reveal counterparts of the Hi- 
cmitting shells. In both cases, the extent of the FIR shell 
is comparable to that of the H I shell, but the two trac- 
ers show notable differences, perhaps due to effects such 
as ionization of the H I from shocks or local destruction 
of dust grains. Both Y UMa and Y CVn also exhibit 
"notch" -like depressions in their Hi and FIR intensity 
along the leading edge of their H I shells and running 
parallel to their direction of space motion. These notches 
may be signatures of instabilities along the leading edge 
of the CSE. While the Hi shells of Y UMa and Y CVn 
are similar in terms of their sizes, morphologies, and line 
widths, the Y UMa shell contains only half the Hi mass 
as that of Y CVn and also shows a distinct rotationally 
symmetric pattern in its H I velocity field whose origin is 
unclear. 

The two highest velocity stars in our study (RX Lep 
and TX Psc) were found to have extended gaseous wakes 
trailing the motions of the stars through the Galaxy. 
Such wakes are formed by the effects of ram pressure on 
the circumstellar ejecta as the stars move supersonically 
through the ISM. The wake of RX Lep (which stretches 
^0.36 pc) is more highly collimated than any Hi wake 
detected to date, a fact that is likely due to a combi- 
nation of the relatively early evolutionary status of the 
star on the TP- AGB and the high density of its inter- 
stellar surroundings. For TX Psc, the orientation of the 
wake (which stretches ~ 0.6 pc) appears to be affected 
by a local transverse flow in the ISM. This flow could 
also explain a velocity gradient observed across the CSE, 
perpendicular to the direction of space motion. 

We find the oxygen-rich Mira R Peg to be surrounded 
by a peculiar "horseshoe" -shaped H I nebula that extends 
over ~1.7 pc on the sky. This morphology may result in 
part from the effects of local ISM flows. However, the 
CSE morphology close to the stellar position is also pe- 
culiar, with two asymmetric lobes of emission stretching 
roughly parallel to the direction of space motion of the 
star. These lobes differ from the head-tail structures seen 
toward other stars in that the "head" of the nebula lies 
at a projected distance of ~0.1 pc southeast of the star. 
One possible explanation is that the observed H I struc- 
tures are the vestiges of an earlier cycle of mass- loss. 

The two undetected stars in the present sample 
(IK Tau and AFGL 3099) are both relatively cool Mira 
variables with high mass-loss rates [M « (4.2, 13) x 
10 _6 M Q yr" 1 ]. Both stars were also undetected in previ- 
ous single-dish studies. For these stars, hydrogen in the 
circumstellar environment may be too cold to detect via 
21-cm line emission or may be largely in molecular form. 
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APPENDIX 

AGE ESTIMATES FOR THE DETACHED Hi SHELLS BASED ON SIMPLE NUMERICAL MODELS 

As a mass-losing star becomes progressively surrounded by an expanding shell of circumstellar material, the collision 
of the supersonic stellar wind with the surrounding medium (local ISM and/or debris from earlier episode of mass 
loss) produces a slowly expanding shell of denser material (Lamers & Cassinelli 1999). Young et al. (1993b) developed 
a model to simulate the formation of these "detached shells" as seen in infrared emission, and Libert et al. (2007) 
subsequently adapted this type of modeling to circumstellar hydrogen. 

In this picture, a detached shell results from a stellar outflow that is abruptly slowed down at a termination shock. 
This first shock defines the inner boundary of the detached shell (ri). The outer limit (7-2) is defined by the leading 
edge (bow shock) where external matter is compressed by the expanding shell (see Figure 5 of Libert et al. 2007). 
Between these two limits, the detached shell is composed of compressed circumstellar and interstellar matter separated 
by a contact discontinuity at a radius r/. The circumstellar matter is decelerated and heated when crossing the shock 
at 7* 1 . As it cools from a temperature (at ri) to Tf (at tv), the expansion velocity decreases, from to Vt, and the 
density increases from to nj. Adopting an arbitrary temperature profile, one can numerically solve the equation 
of motion between r± and r / . 

Here we have applied a similar method to modeling the H 1 shells of RX Lep, Y UMa, Y CVn, and V1942 Sgr. 
Models for three of the stars (RX Lep, Y CVn, and V1942 Sgr) have previously been published based on single-dish 
Hi observations (Libert et al. 2007, 2008, 2010a), but our new imaging observations provided improved constraints 
on the model parameters. 

We have not attempted to compute models for either TX Psc or R Peg. In the case of TX Psc, the relative high space 
velocity and large fraction of the circumstellar ejecta present in its extended trailing wake makes it a poor candidate 
for a stationary and spherically symmetric model. For R Peg, the complex H 1 morphology and line shape also imply 
that a simple model fit of the type described here is unlikely to be meaningful. 

Following Libert et al. (2008), the temperature profile in the detached shell of the four modeled stars is constrained 
by the observed Hi line profile. For the external part of the detached shell (between 77 and r 2 ), we assume a density 
profile that falls off as 1/r 2 , with the condition that the total mass (Mdt.cs) is given by the equivalent ISM mass 
enclosed in a sphere of radius ra< Inside ri, we assume a spherical wind in uniform expansion, whose properties 
(mass- loss rate, outflow velocity, and systemic velocity) match those derived from CO observations (Table 1). We 
assume mass lost from the star comprises 90% H atoms and 10% He atoms by number (Glassgold & Huggins 1983). 
The outer radius of the shell, r*/, is constrained by the VLA images. 

Based on our models, the detached shell is expected to produce a narrow, quasi-Gaussian emission peak at or near 
the stellar systemic velocity in the spatially integrated spectrum, while the freely expanding wind should produce a 
weaker, broader line component with peaks located on either side of the stellar systemic velocity, near ± V a (see 
Le Bertre et al. 2012). In general, the shell component dominates the global spectral profile. 

Our best- fit parameters for the four stars that we have modeled are summarized in Table Al, and model fits are 
overplotted on the observed spectra in Figure In all cases, our toy model provides a reasonable fit to the strength 
and width of the global H 1 line profile. 

For the present study, the most useful outcome of our models is an estimate for the age of each of the shells, which 
can be obtained from a sum of the parameters t\ and t;ds in Table Al. Not surprisingly, we find the smallest age 
(^90,000 yr) for RX Lep, a star that is believed to be in the early stages of the TP-AGB (see §4.2[) . For the other 
three stars, our models imply that the Hi data are tracing more than 10 5 years of stellar mass loss history. In the 
case of Y UMa, our derived shell age is comparable to the age of the FIR-emitting shell previously derived by Young 
et al. 1993b (9.2 x 10 4 yr), while in the cases of the carbon stars Y CVn and V1942 Sgr, our new values are 6.8 and 
6.6 times larger, respectively. 

To obtain a more complete picture of these four stars, future models will ultimately need to take into account the 
motion of the stars through the ISM and the resulting distortions that this produces on the Hi shells. However, this 
will require sophisticated hydrodynamic simulations that are beyond the scope of the present paper. Furthermore, 
we note that our model fits break down at small radii (r <ri), systematically predicting more flux than is observed 
in both the shell component (central peak) of the line profile and the freely expanding wind ("horns"). This may be 
due to an uncertainty in the temperature profile adopted in our models, short-term fluctuations in the mass-loss rate, 
asymmetric or non-spherical mass loss, and/or a breakdown in our assumption that the wind is largely atomic as it 
leaves the star. These possibilities will be investigated in a future work. Fortunately the ages that we derive for the 
detached shells (Table Al) depend mostly on their large-scale properties and are relatively insensitive to these details. 
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Fig. 28. — Global (spatially integrated) Hi spectra for four stars from the current sample (solid histograms), with spectra derived from 
numerical models overplotted (dot-dash lines; see Appendix A): (a) RX Lep; (b) Y UMa; (c) Y CVn; (d) V1942 Sgr. The model spectra 
have been smoothed to match the spectral resolution of the observations. For RX Lep, the spectral peak blucward of the stellar emission 
results from interstellar contamination. 



TABLE 1 

Detached Shell Model Parameters 



Parameter 


RX Lep 


Y UMa 


Y CVn 


V1942 Sgr 


M (Mq yr- 1 ) 


2.0xl0~ 7 


2.6xl0~ 7 


1.7 x 10~ 7 


1.0 x 10~ 7 


A 1 


1.3 


1.3 


1.3 


1.3 


*l (yr) 


2.0xl0 4 


3.6 xlO 4 


2.8xl0 4 


2.4 x 10 4 


*ds (yr) 


7.0xl0 4 


8.8xl0 4 


4.2xl0 5 


7.0 x 10 5 


n (pc) 


0.087 (2.'0) 
0.11 (2/6) 
0.17 (4') 


0.22 (2/0) 
0.26 (2/3) 
0.34 (3/0) 


0.22 (2/8) 
0.32 (4/0) 
0.40 (5/1) 


0.12 (0/8) 
0.22 (1/4) 
0.44 (2/8) 


Tf (pc) 


r 2 (pc) 


ToCsTf 1 ),^ (K) 


20, 528 


20, 1071 


20, 1804 


20, 746 


T f (=T 2 )(K) 


157 


221 


173 


87 


vo(= v^),vf (km s 1 ) 


4.2, 1.07 


6.0, 1.52 


7.8, 1.96 


5.0, 1.27 


Vf (km s _1 ) 


0.16 


0.18 


0.07 


0.04 


V2 (km s _1 ) 


1.35 


1.2 


1.0 


1.26 


jri^ (H cm -3 ) 


14.0, 54.8 


1.9, 7.55 


1.0, 3.8 


2.8, 11.2 


nj 1 ,n~j (H cm~ 3 ) 
ri2 (H cm -3 ) 


228, 2.3 


45.0, 0.5 


50.0, 0.6 


119, 1.6 


1.0 


0.2 


0.4 


0.4 


M r<ri (Mq) 


4.0 x 10~ 3 


9.5 x 10~ 3 


4.8 x 10~ 3 


2.4 x 10~ 3 


M VT ,CS (Mq) 


1.4 x 10~ 2 


2.3 x 10~ 2 


7.17 x 10~ 2 


7.0 x 10~ 2 


M D T,EX (Mq) 


7.7 x 10~ 4 


1.1 x 10" 3 


2.4 x 10~ 3 


7.3 x 10~ 3 



Note. — Notations follow those adopted by Libert ct al. 2007. Matter lost is 
assumed to be 90% neutral hydrogen, 10% helium by number density. Distances are 
adopted from Table 1. 
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